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A. BACKGROUND 


Since the introduction of the sampling theorem and tne 
matched fiiter, digital communications techniques have 
developed into a highly proficient discipline. The marriage 
of this discipline with the rapidly expanding space fprogran 
has cesulted in ccmmunication sateilites employing a 
Multitude of digital signal modulation tecaonigues. A 
modulation technique is a method of transmitting the 
information contained in a message by varying cr modulating 
the characteristics of a carrier waveforn. These metnods 
offer a range of advantages and disadvantages depending on 
the specific characteristics of the modulation technique 
employed. The applications of the various Signal modulation 
tecanigues likewise vary. in’ order to understand this 
exciting new field it is necessary to look at some of the 
aspects of satellite Communications systems in general. 
Then an investigation will be made into the general 
attriktutes of digital communications and their relationship 
to the satellite systen. . 

The abiilaty £0 understand these digital Signal 
modulation techniques is the first step in being able to 
intercept, identify and demodulate unknown digital signals. 
These digital signals, transmitted from unknown sources, are 
believed to display frequency characteristics peculiar to 
the mcdulation technique employed in the encoding process. 
The digital computer offers unigue opportunities in 
Simulating these modulation techniques, in siynal processing 


and in decoding of an intercepted signal. 


AA 


B. SPECIFIC GOALS 


It is the specific goal of this thesis to investigate 
the open literature on the topic wo digital sig@ad 
modulation technigues in satellite communications. This 
includes a basic understanding of the communications 
satellite system as well as the specific techniques employed 
in signal modulaticn. Additionally, once a basic 
understanding of these digital signal modulation technigues 
is achieved, computer code in the FORTRAN programming 
language will be developed which simulates these modulated 


Signals. The statistics of the time-varying Fast Fourier 
Transforms (FFT) of these Simulated signals will be 
lavestigated. The purpose of this analysis wiil be to iead 


to fcllow-on research in the area of Signal analysis’ of 
intercepted digital signals whereby they can be classified 
by their FFT as employing a specific .digital signage 
modulation technigue. 


C. SCOPE OF THE PROJECT 


Althougon an indepth analysis of all the digital Signal’ 
modulation technigues which will be introduced involves 
considerable higher level mathematics, it is not-within the 
scope ot this project to delve heavily into the mathematics. 
It would be advantageous, however, for the reader to have 
had integral and differential calculus and an introduction 
to statistics. Also a course in electrical engineering nay 
prove helpful but is not essential since in actuality it is 
the electrical engineering aspects of satellite 
communications that this tutorial is attempting to present. 
The computer programming will be accomplished utilizing the 
techniques of software engineering and top down modular 
design. Existing blocks of code will be used as modules 


whenever appropriate and available. Again it is Gem 


ultimate purpose of this project to examine a variety of 
digital signal modulation techniques and develop computer 
code that simulates common Signals used in satellite 
communications that have been produced by one of the many 
digital signal modulation techniques to be investigated. 
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IIT. UNDERSTANDING SATELLITE COMMUNICATIONS 


A. HISTORY AND APPLICATIONS 


Forty years ago, in 1945, Arthur Clarke first envisioned 
the use of space stations placed in geosynchronous orbit for 
communicating to different points on the earth [Ref. 1]. 
Nine years later in 1954, J.R.~ Pierce of Bell Lakoratories 
performed a system analysis on such a communications system 
[Ref. 2]. In 1957, the launch of Sputnik demonstrated the 
feasibility of using a satellite for’? Fast such an 
application. However, by 1961 the only satellite 
communications tecknologies which had been demonstrated were 
the Ccurier 1B satellite, a short-life active retransmission 
teletype communications satellite ina 1000 km orbit, and 
the Echo I balloon ina 1600 km orbit which demonstrated 
passive reflection of powerful microwave signals from one 
earth station to another. Active microwave communications 
demonstrated in Projects Telstar and Relay were years away 
fRef. 3]. The first geostationary orbit was achieved by 
Project Syncom in 1963 [Ref. 4]. 

In 1964, communication organizations from several 
countries joined together to form the international 
organization of INTELSAT (International Telecommunications 
Satellite Organization). INTELSAT"sS purpose was to develop 
a satellite network which would provide truly glottal 
communications capabilities. This resulted in the launch in 
1965 of the world's first commercial communications 
satellite, INTELSAT I, also known as “Early Bird". With 
"Early Bird", telecommunications utilizing satellite relay 
were established between the United States and Europe. 
{Ref. 3] 
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The reliability of these satellite systems has improved 
dramatically since INTELSAT I in 1965. Reliability of 
individual links in the system approach 99.99 percent or 
higher. Total system reliability exceeding 99.9 percent is 
common [Ref. 5], making satellite communications more 
reliable than many other modes of communications. in thas 
sense, reliability is a measure of the probability that no 
failure will occur in a respective channel or in the systen 
during the design life of the -satellite. 

In crder for the operability, capability and reliability 
of these systems to have developed at this pace, it was 
necessary for the technologies associated with them to 
develop as rapidly or even more rapidly than the systems 
themselves. Engineering sciences and specifically the 
fields of aerospace and electrical engineering historically 
have reguired between 7 and 10 years to take a concept from 
operational requirement to full scale operation. This was 
not the case with the concept of communications satellites 
which has seen four generations ot INTELSAT sateilites 
within a decade's time. [Ref. 3] 

The “Early Bird" satellite weighed approximately 38 
kilograms and possessed limited power and bandwidth capacity 
enabling LE to carry only 240 two-way telephone 
conversations. Today's communications satellites represent 
order-of-magnitude improvements in many important oferating 
parameters such as fower and bandwidth. Addrezonally, 
increased effective radiated power from the technigues of 
Stabdilized earth pointing antennas have greatly increased 
the capacity of later generation communications satellites. 
EVTELSAT V, the current generation of communications 
satellites, is a three-axis stabilized platform using not 
only the 6/4 GHz frequency band (6 GHz receive/4+ GHz 
transmit) as in earlier generation satellites but also a 500 
BHz kandwidth available in the 14/11 GHz band. The 


UES) 


separation in the receive and transmit frequencies is 
necessary to prevent interference during simultaneous 
operation ot the receiver and transmitter. Ancther factor 
increasing the capacity of present generation satellites is 
the increase in primary power available in the satellite. 
{Ref. 3] 

The technologies which have contributed to the evolution 
of communications satellites come from two primary sources, 
hamely technology from the Space program of the 1960's 
supported by NASA and DoD and communications technology due 
largely to commercial and private sector contributions. 
{ Ref. 3] 

Electronic devices and components have contributed 
Significantly to the rapid development of satellite systems. 
These electronic devices range from something which is now 
considered basic, i.e., the transistor, to devices such as 
Traveling Wave Tube amplifiers, lightweight antennas and 
antenna feed systems. There have been major improvements in 
satellite power sources including more efficient solar ceils 
and high storage Capacity/lightweight batteries. 
Additionally, two Significant aevelopments in electrical 
engineering have made modern day digital signal processing a 
reality. They are the sampling theorem and the matched 
filter, {Ret. 3] 

Communications satellites of the future are likely to 
utilize onboard signal processing. Sigtal processing 
functions such as stgqnal reshaping, Switching and/or 
multiplexing and compression could soon take place on the 
satellite due to the reduction in size and weight of the 
necessary hardware. Operation at over 100 megabits per 
second are envisioned. Onboard signal processing will 
greatly reduce the expense and complexity of earth stations 
thereby making services of a satellite available to a wider 


range of small and gecgraphically disperse users. [Ref. 3] 
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Satellite communications hold great vromise to provide 
service to a multitude of users over a wide area. 
Applications lie not solely in retransmission but also in 
data collection frcem that same large geographic area. 
Military users have definite applications in these areas 
when warning, intelligence and surveillance systems provide 
digital inputs to a central source. Although most of the 
present applications for satellite communications are still 
provided by the government, commercial applications nave 
seen tremendous increases in the last several years. 
{Ref. 3] 


Be ORBITS AND LIMITATIONS 


Satellite orbits are generaily categorized as either 
equatcrial (0 degrees inclination), polar (90 degrees 
inclination) or inclined at some angle other than 0O or 90 
degrees relative to Lhe “Spin “axis ”* Of the earth as 
morustrated in Figure: 2.1 (Ref. 6]. Each satellite orbit 
has a characteristic velocity which is dependent on the 
height of the orbit and the orbit's eccentricity. 
Eccentricity iS a measure of the degree to which the orbit 
approximates a circle. A circle has eccentricity equal to 
zero since the focus is located at the center. See equation 
2.%*and figure’ 2. 2; 

A communications satellite in elliptical orbit about tne 
earth obeys Kepler's second law of planetary motion. That 
is, a satellite's constant angular momentum about the earth 
Means that its areal velocity also remains conrstart. See 
Prairies ~ 2. 36 Of particular interest is the satellite 
velocity at apogee (Va) and perigee (Vp) given by equations 
wa22 ana~ 2.3. 

Note that for a circle, e = 0 and Ra = Rp. Therefore a 
Sateigate fn “circular “orbit about the Garth ~has°an* orbital 


velocity as given in eguation 2.4. 
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Figure 2.1 Satellite Orbits 


Most communications satellites are placed in a circular 
equatorial orbit where the desire is to stabilize them over 
a fixed point on the surface of the earth called the 
subsatellite point. This type ef. .Orbi £ is... called 
geostationary or stationary. Any satellite which has an 
orbital period equal to the period of rotation of the earth 
is called synchronous or geosynchronous. These terms are 
used almost interchangeably in most literature. As 
mentioned for a geosynchronous communications satellite, the 
orbital period of the satellite, TT, must be equal to the 
period. of. otatLon of. the. earn. Kepler's third laweoe 
planetary motion can be rewritten to yield equation 2.5. 

The period of revolution of the earth for a sidereal day 
is 23 hr 56 min 4 sec. For that given period there is only 


one satellite altitude as expressed by Kepler's third law. 
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Figure 2.2 Eccentricity 


eccentricity (e) = c/a (eqn 2.1) 


distance from focus 


= semi-major axis 


b = semi-minor axis 
Ra = radius cf apogee 
Rp = radius cf perigee 


By rearranging terms that altitude is given in equation 2.6. 


Simce the, erbhi-t.1s Girculatr,,-. ai=. Ra.= Rp,and-the height of 
the orbit above the surface of the earth is o = a —- Re; or 
35,804 kao. 


One disadvantage of a geosynchronous communications 
satellite is the lack of global coverage. These satellites 


provide excelient coverage of the subsatellite point and 
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igure 2.3 Kepler's Second Law of Planetary Motion 


(k/Ra (1-e) ) **.5 


(k/Rp (18) 4"s5 
GMe = 3.98866 x 10**-11 m3/s2 
universal gravitational constant 
6.67 x 10**-11 N m2/kg2 
newtonrs = mn kg/s2 
mass of the earth = 5.98 x 10**24 %kg 


(k/Re) =*25 


= Re + hb 
radius of the earth = 6.37 x 10**6 
satellite orbital altitude 


2 pi at*ic5 7 k=e25 


period of rotation 


semi-major axis 


(eqn 2.2) 


(eqn 2.3) 


{eqn 2.4) 


{eqn 2.5) 


See Cle 2 pa) Stes 3S (k) FF 173 (eqn 2.6) 
a = 42,173 ko 
laterally to a latitude of about +80° {Ref. 6]. This is 


generally no problem for commercial applications since the 
polar regions do not require a significant degree of access. 
The military, on the other hand, does have an interest in 
communications in the polar region and therefore has several 
Communications satellites that have orbits inclined at 
various angles. This type Of OEbit generally has 
disadvantages of lack of continuous coverage and a much more 
complicated systen of ground tracking and receiving 


staticns. 


C. FREQUENCY BAND CONSIDERATIONS 


Although there appears to be an infinite nuwber of 
frequencies available for communications, restrictions do 
exist as to those which are practicable. Limitations on 
available frequency bands for satellite communications ere 
due tc the need to select segments of the electromagnetic 
spectrum which reduce noise and interference and are most 
favorable in terms of power efficiency and propagation 
distortions. Trade offs must be made to arrive at the 
optimum frequency for a particular application since single 
frequencies seldom offer the best performance for ali 
variables. The problems arise when consideration is given 
to the number of users requiring the same frequency bands 
including terrestrial communications networks. Since the 
problea of interference is glotal, a worldwide organization 
has been established to assign frequency bands for various 
applications. This organization is called WARC, World 


Administrative Kadio Conference [Ref. 6]. Illustrations of 


aA 


the pcrtions of the electromagnetic spectrum in guestion and 
current allocations of satellite frequency bands ate shown 
in Figures 2.4 and 2.5 [Ref. 6]. 
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Figure 2.4 Electromagnetic Spectrun 


In general, as long as the hardware and technology will 
support it, higher frequencies are more desirable and more 
in demand because they offer higher theoretical capacity. 
This is due to the fact. that only a percentage cf sytae 
Carrier frequency is capable or actually transmitting a 
Signal of a given bandwidth. Higher frequencies wouid also 
experience less interference with existing land and 
satellite systems. A further discussion of the bandwidths 


available as a function of frequency will be included at a 
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Allocated satellite frequency bands, WARC 1979 (7 = uplink, 
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Figure 2.5 Satellite Frequencies 


aber) pornt qeunders) thes topic» of actual.« digital sagnat 
techniques. 

As an additional note, in order to assist in the ability 
to handle nore Signals in an obvicusiy limited 
electromagnetic spectrun, several muitiple access schetes 
have reen developed. The mest common tecnnigues involve 
rrequency-division multiple access (FDMA) where the 
allocated satellite frequency band is divided among the 
users into specific uplink and downlink frequencies. Next 
is time-division multifle access (TDMA) where the sateliite 


frequency band is shared by all users by carefully dividing 
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the time any one user has access. Finally there is 
code-—division multiple access (CDMA) which involves 
modulating a specific address code which has been 
superimposed on the signal directly onto the carrier. In 
this manner all users share the satellite frequency band and 
only those receiving stations that can demodulate the 
address code can receive the specific signal. [Ref. 6] 


D. WHY DIGITAL MODULATION 


In general, digital signal modulation techniques are 
superior to analog signal modulation technigues used in 
satellite communications for the following reasons: 

1. Compatibility with digital computers 
2. Economic advantages 

3. High degree of flexibility 

4. Less susceptible to interference 


5. Quality of signal independent of transmission 
distance and network makeup 


Clearly one of the most distinctive advantages of 
digital technigues over analog techniques involves computer 
‘ applications, Properly formatted digital signals can be 
used to represent any analog signal (more on this under the 
discussion of the sampling theorem). Once in the digital 
format these signals can be easily manipulated within the 
digital computer. Arithmetic operations can Le applied as 
well as logical operations. The signal can te _ stored 
without alteration or delayed and theretore can be used to 
"simulate" real physical situations. 

The hardware associated with digital circuits is 
free from drift or aging and does not require calibration. 
Additionally digital circuitry is compatible with present 


day integrated circuit technology allowing a standardized 
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building block construction approach. The operating 
characteristics of these systems employing ae adaqa tal 
computer can be changed by altering the software rather than 
the hardware as is the case fcr analog systems. Finally, 
the use of digital computer technology allows time 


multiplexing. 
2. Flexibility and Economy 


The flexibility and economy of digital satellite 
communications comes from the fact that more and nore 
processing can be done onboard. This allows the uplink and 
downlink to be compietely separated. This regenerative 
Mature makes it possible for low error rates and aigh 
Peladbiiity through. the use of digital techniques not 
available to analog systems. Because digital signal 


processing or multiplexing is less costly than for analog 


Signals, Simpler and cheaper interfaces between earth 
staticns and terrestrial communications networks are 
possitle. Additionally, there are reduced production costs 


and increased capacity associated with digital circuits. 
(Ref. 7] 


3. Quality and Interference 


The capability of digital systems to regenerate tie 
Signal and allow for multiple switching and Signai 
processing without degradation in signal quality makes the 
digital Signal basically independent of transmission 
distance. Multiple hops from satellite to earth staticn or 
Satellite to satellite are possible without accumulation of 
the ncise characteristic of analog systems. Additionally, 
digital systems are capable of operating at a Signal to 
noise ratio of 20 dB to 30 dB as compared to analog systems 
requiring a much more powerful signal. [Ref. 8] 
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Essential to the understanding of digital | signad 
modulation techniques are a few basic tools of the 
electrical engineer and the mathematician. These tools will 
be developed and elaborated on to the extent necessary to 
understand their applicability to the subject of digital 
communications. The description is not meant to be a 
detailed investigation of the respective topics. where 
relevant, the application of the concept being described 


wili te mentioned. 


A. FCURIER TRANSFOR SH 


In mathematical terms, voltages can be expressed as 
functions of time or of frequency. It is more commcn to see 
voltages represented as a function of time as in equation 
J eelie 


v(t) = A cos(wt) (eqn 3.1) 


> 
iH} 


amplitude 


angular freguency = 2 pif 


frequency = 1/T 


period 


It 1S important tc note that both the time and frequency 
functions are representations of voltage and as such may be 
used interchangeably. The Fourier representation Cv (eh f= 


V(f) ais avoltage descriptor in the frequency domain {a 


function of frequency) while v(t) is a voltage descriptor in 
the time domain (a function of time). They are different 
but eguivalent and either voltage descriptor may be used, 
depending on the concept being explored, to best represent 
the voltage within context. 

The Fourier transform is of principle interest rather 
than the Fourier series since the latter is applicable only 
to periodic voltages. The Fourier transform is applicable 
to voltage pulses, random voltages and other rnon-periodic 
voltages. 

DEFINITION: 


CO 
Frv(t)] = Vw) = [v(t) ex*-jut dt (eqn 3.2) 
. -o 
V(t) <-------------> V(¥) 
time domain frequency domain 


Remembering Euler's formula 
e**-jwt = cos(wt) - j Ssin(wt) (eqn 3.3) 


the Fourier transform becomes 


Co 
At v(t) ] = Vw) = f v(t) [cos (wt) - j sin(wt)]}© (eqn 3.8) 
fore) woe CO 
Vv (w) = [ v(t) cos (wt) - 3 [ v(t)sin (wt) (eqn 3.5) 
—oo — oo 


Since it will be seen that digital communications deals 
Pplrigarily with pulses of finite duration (expressed as 
period, TD) it 1S worthwhile to examine the Fourier 
transform or a pulse cf amplitude A and duration T. 

eK 2h EK Tee 
let v(t) = 


0; elsewhere 
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Figure 3.1 Square Pulse 


Figure 3.1 iS a representation of v(t) or the voltage 

expressed in the time domain, The position of v(t) on the 
4 

t-axis was chosen for convenience of integration but could 


have been Situated anywhere on the time line. 


Fiv(t))] = vi) = felt) ete—gue at : (eqn 3.6) 
-% = Ts. 4 ; 
vu) = Jere jut dt + {A e**-Jjwt dt + (eqn 3.7) 
0se**-jwt dt ve 
T/2 Tyan 
V(w) = fa e**-jut dt 2 (eqn 3.8) 
te 


The integration above may te attacked either nead on or 
by substituting cos(wt) - j sin(swt) for e**-jwt. The diver 
approach is illustrated due to the relative simplicity ot 
the integrand. 


y 


z 
V (w) =fa e**- jwt dt (egn 3.9) 
Li 
122 
V(w) = -A/jwle**-jwt] | (eqn 3.10) 
-1T/2 
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V (w) -A/jwle**-jwT/2 - e** jwT/2) (egn 3.11) 


V (w) A/jwl e** jwI/2 - e**-jwT/2 ] {eqn 3.12) 


By substituting 2pi f = w, the following result is obtained: 


¥V(£) = A/j2 pi fL[e**j2 pi f£T72 - (egn 3.13) 
e**-j2 pi f£T/2] 


V(£) = A/j2 pi fl e**j pi £1 - e**-5 pi FT] 


Using Euler's formula, ise., sin 9 = (e**j8-e**-j9) /2j: 
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V(f) = 2jA/j2 pi f£[ (e**j pi fT - “(eqn 3. 14) 
e*F— 7) pi £1) /2 7) 
MeL = A/pt FE Sint pi fT) 7] 

Knowing that sin x / x = sinc x 
Tce) = AL/pi- LT Sin (pi £1) | (egn 3.15) 
ec) = AT Sipe (pi tT) 


THe Si nGwluncrion..1s common ip.digital,. electronics. and 
prots.as. thes product, of. .sin( pi. iT) and 1/4p1 fT) aS aan 
Figure 3.2. In Figure 3.2, V(f£) in the frequency domain is 
equivalent to v(t) in the time domain. Note that as the 
pulse]. 14, Gets longer, 1/1 gets smaller or the first zero 
crossing of the sinc function occurs ata lower and lower 


frequency. 
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Figure 3.2 Plot of Sinc Function 


1. Amplitude and Phase Spectrum 


Recall that V(w) can be represented as in equation 
Je. Oc The cosine term is the real part while the sine term 
is the imaginary part. By referencing Figure 3.3, a brief 
review of the complex ‘plane is accomplished and its 
relationship to the amplitude and phase spectrum of a given 


voltage is'represented. See eguations 3.17 through 3.20. 


Co (oe) 
V(w) = [v(t) cos (wt) = ify (t) sin (wt) (eqn 3.16) 
=e real — COimaginary 


{V(f) 1 as called the amplitude spectrum or the given 
voltage. The amplitude spectrum can also be calculated by 
using the complex conjugate or the Fourier transform andiis 
always positive aS shown in equation 3.21. The phase 
spectrum, 986, of the function in question is represented by 


the arctan{ imaginary/real] and is illustrated in Figure 3.5. 
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Figure 3.3 Phasor Diagran 


r= {[a2 + b2] (eqn 3. 17) 
ao * jbe=—Fa- + bhej*+.5 e** 756 (eqn 3.18) 
v(£) = {real?@ + imaginary? ]**.5 e**j6 (eqn.3., 19) 
AV €£) | = [real2 + imaginary? ]**.5 (egn 3 203 
since |e**j@| = 1 

Perle = 1 ¥ (fyevs (2) jF*. 5 (eqn 3.21) 


Several properties of the Fourier transform are 
useful in the Study Of digital signals. They are 
represented here without proof and without a great deal of 
detail. 

a. Linearity 
PEW (ey) <<>> Vil) and 
af Vert) <-> V2(f) then 
Flavi(t) + bv2(t)] = avi(f) + bv2(f) 
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Figure 3.4 Amplitude Spectrum of Square Wave 
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Figure 3.5 Phase Spectrum of Square Wave 


b. Time-delay 
Giiv (t-to) J = V(f) e**-jwto 
Note that the amplitude spectrum of the delayed 
verSicn is the same as the amplitude spectrum of the 
undelayed version. It is comforting to note that aaa. 
Fourier transform of a pulse is the Same tomorrow as it is 
today. 
c. Scale change 
Ftv (at) ] = What v(£/a) 
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d. Frequency translation 
v(t)cos(2 pi fe t) <-> 1/2[V(f+fc) + V(f-fc) ] 
v(t) is any voltage 
cos(2 pi fc t) is a carrier wave of frequency fc 
Since understanding of this very important property 
of the Fourier transform is essential to the understanding 
of digital signal modulation it is expanded slightiy here. 
As we have seen, the Fourier representation of a 
voltage pulse of amplitude A and duration T is the sinc 


function of amplitude AT as illustrated in Figure 3.6. 
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Figure 3.6 Square Wave in Time/Sinc Function in Frequency 


The translation is the product of v(t) and in this case 
SOst2 pl cc tj.” In the time domain this product only exists 
in the interval between -T/2 and T/2 as in Figure 3.7. a 
Significance of this translation and its relationship to the 
bandwidth of tne voltage will be discussed further in the 
section dealing with bandwidth. 
e. Differentiation 
d v(t) /dt <-> jw V(f) 
A differentiator could be used as a clock for timing 


but would never be used in the presence of noise since the 
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Figure 3.7 Multiplication and Translation Diagrams 


result is exaggerated in the presence of high freguencies 
because w= 2 pi f. 
f. Integration 
v(t) dt <-> 1/w V(£) 
An integrator could be used to reduce the effects of 
noise since w=*2 pi f is in the denominator tending to 


deemphasize the presence of high frequency noise. 


Be. THE SAMPLING THEOREM 


Essential to the uaderstanding of digital communications 
is the sampling theorem which was first introduced by 
Nyquist in 1928 [Ref. 9], and later by Shannon in leeee 
{Ref. 10]. The sampling theorem states that any voitags can 
be uniguely represented by appropriately spaced sampie 
values of the original voltage. More correctly stated, the 
Sampling tneorem places limits on tne accuracy with wnich a 
Signal can be represented. 

The implication is that an analog signal can be 
represented digitally or by a set of numbers, i.se., Samples. 


A description of the sampling theorem follows. 
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Given any analog voltage, v(t): asin Fig ures3'23,- the 
sampling theorem says that the entire analog signal is not 
required to accurately represent the voltage but only 
samples of it, call them vs(t). Figure 3.9 shows samples of 
a representative analog voltage. Samples can be taxen of 
v(t) at every T seconds for a period of seconds. This can 
be accomplished by the use of a voltage clock, call it 
we(t). The sampling can be viewed graphically as a block 
diagram representing an analog voltage multiplier as shown 
Py Sagure. 3.10. To be of further use in the understanding 
of digital communications we are interested in a frequency 
description of the sample voltage, vs(t). Note that the 
system which describes the obtaining of vs(t) involves a 
voltage multiplication. It was demonstrated in the 
proceeding section that voltage multiplication amounted to 


frequency translation, a property of the Fourier transforn. 
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Figure 3.8 Representative Analog Voltage 


First it is necessary to find the frequency description 
of the clock, vc(t), Let vc(t) be a periodic square wave of 
height 1 and duration das in Figure 3.11. Sance vco(t) 2s 
periodic, ew Can be shown ana the Fourier series 


representing vc(t) is given by equation 3.22. 
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Figure 3.9 Sasples of a Representative Voltage 


v(t} ——>+} AVM | vs (+) 
tp | 


ve (+) 


a 


Figure 3.10 Analog Voltage Multiplier 


Again it is noted that vs(t), the sample voltage, is the 


product of v(t), the original analog voltage, times vc(t), 


the clock voltage. In other words vs(t) = v(t) times a 
series cf cosine terms. 
(o) 
vc(t) = ao +) an COS(2. pi By te. th (eqn 3.22) 
ns 


ao and an are left unevaluated 


, 


ve (4) 
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Figure 3.11 A Representative Voltage Clock 


(ore) 
vs(t) = v(t) [ao +.) an cosi(Z pi on fc t) ] (eqn 3.23) 
h=1 


Now, assuming any form of the Fourier transform of V(t); 
as in Figure 3.12 and since vs(t) isa “product and. by 
utilizing the frequency translation property of the Fourier 
transform, the following is obtained as a representation of 
the frequency spectrum of vs(t). See Figure 3.13. 

Remember that the curve highlighted in the box in Figure 
ao13-ts-the—rourier represention- “of ~-v (tt), the original 
Signal. This original signal can now be recovered by 
filtering with an appropriate low pass filter of bandwidth 
equal to or greater than B. This low pass filter wceculd only 
permit the reception of that portion of the signal which 
represents the original voltage. 

The only question left to resolve is how often to take a 
sample. Again referring to the diagram in Figure 3.13, it 
is noted that in order to prevent any overlap of successive 


translations (aliasing) 
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Figure 3.12 Fourier Transform of v(t) 
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Figure 3.13 Fourier Transform of vs (t) 


fOr = 5 >> a8 

£0) D>. 258 
or the frequency of the clock, fc (the sampling rate = 1/T) 
must be strictly greater than 2 times the largest 
Significant Fourier frequency component present. Since B 
will vary with different v(t), the sampling rate necessary 


to uniquely recover that v(t) also varies. 
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In summary on the sampling theorem, it can be said that 
any analog signal can be represented as certain sample 
values spaced the appropriate distance apart. The sampling 
theoren places limits on the accuracy of this 
representation. These sample values can then be transmitted 
from one position to another using analog to digital 
conversion and any one of a variety of moduiation 
techniques. It has been shown that the Fourier frequency 
representation of the sample is not equivalent to the 
Fourier frequency representation of the source voltage; 
however, the source voltage can be recovered at tae 
receiving erd by filtering. A block representation of the 


system is shown in Figure 3.14. 
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Figure 3.14 Analog Voltage Multiplier and Low Pass Filter 


C. THE AUTOCORRELATICN FUNCTION 


The above concept certainly seems simple enough. ik. a 
Signal is present and it is desired to transmit it froa 
point A to point B, all that has to be done is to take 


appropriately spaced samples of the signal, encode then, 


a8 


somehow transmit them to point B, filter the received signal 
and the unigueness of the original voltage has been 
reproduced. Unfortunately it isn't quite that easy. The 
reason it is not, is due to the presence of noise. The 
sources of noise will not be discussed here, however, the 
effects of noise in general terms and how a faint signal can 
be recovered in the presence of noise will be discussed. 

The time varying descriptor of voltage and the frequency 
varying descriptor of voltage have been introduced. Both 
are precise mathematical descriptions of something (voltage) 
that is deterministic. That 25, it can be described in 
mathematical or graphical terms during any period of time. 
Such is not the case for noise corrupted voltages which are 
entirely randon. Therefore other descriptors of voltages 
must be employed in the presence of noise. They are partial 
descriptors of voltage since a random signal cannot be 
described with precision. These partial descriptors are the 

1. Autocorrelation function and the 
2. Probability density function (p.d.f.) 

It iS important here to note the difference between the 
source voltage at the receiver and the sample voltage. The 
source voltage at the receiver, vsr(t), is the digitally 
converted sample voltage, vs(t), modulated onto a carrier by 
one of the digital modulation techniques yet to be 
discussed. Understanding the autocorrelation function is 
the basis for understanding how a decision is made that 
vst (t) is present at the receiver in the presence of noise. 
It should be remembered that the exact form of the carrier. 
is kncewn at both the transmitter and receiver. 


First of all, an additive noise model is assumed where 


the signal at the receiver, vr(t), is equal to the signal 
which is being transmitted, vsr(t), (A 5e5 5 the digital 
samples of v(t) modulated onto the carrier), plus random 


noise, vn(t) 
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vr (t) = vsr(t) + vn(t) (eqn 3.24) 


vr(t) = voltage at the receiver 
vsc (t) = Signal voltage at the receiver 
vn(t) = noise voltage at the receiver 


If it is redlized that most receivers operate on the 
principal of detection of DC voltage, the problem becomes 
one of rectification of the received signal and 
determination if the DC voltage, Characteristic of the 
transmitted signal is present. 

A common type of rectification of an analog voltage 
involves squaring the input waveforn. If the given signal 
at the receiver is vr(t) aS in eguation 3.24 above, squaring 
the waveform introduces the square of not only the desired 
Signal but also the square of the noise term. If instead of 
squaring vr(t) and introducing or at least not eliminating 
fEmemnorse, Vr({t) is multiplied by vsr(t), the results in 


equation 3.25 are obtained. 
vo (t)evsr(t) = vsr**2(t) + vsr(t) vn (t) (egn 3.25) 


By averaging, all that remains is vsr**2(t) since the 
average of vsr(t)vn(t) is zero because vn(t) is random and 
the average of any random voltage is zero. The DC component 
of vr(t)evsr(t) ais vsr : Any remaining AC component 
can be removed by a low pass filter. Graphically in block 
diagram form: Chis receiver looks like Figure Ses) 
illustrated below. If vsr**2(t) > 0 then vsr(t) 1S present 
in the signal vr(t). If vr (t) is pure noise or some other 
Signal is present exclusively, then the output of the 


receiver will be 0. 
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Figure 3.15 Rectification through an AVM and LPF 


In practice, however, the exact waveform of vsr(t) may 
not be known or it may be a time delayed or distorted 
versicn of the original Signal when it arrives at the 
receiver. If the distortion or delay is significant enough 
there will be little agreement or "correlation" in the 
receiver, 

The solution is to multiply the received signal by a 
series of time delayed approximations of the transmitted 
signal. The Signal is present whenever the output of the 
series of voltage multipliers or correlators exceeds a 
certain threshold approaching vsr**2(t). The concept is 
illustrated in Figure 3.16. 

For the concept illustrated in Figure 3.16 to work, the 


average value of vr(t)evsr(t-p) must have a DC component. 


vc (t) evsr(t-p) = vsr(t)vsr(t-p) + (egn 3.26) 
vn (t) vsr (t-p) 


vc (t) evsr (t-p) VSI (t) VSE (t-p) + (eqn 3.27) 


vn(t) vsr (t-p) 
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Figure 3.16 Correlation Device 


But the average of vn(t)vsr(t-p) = 0 since vn(t) is random and 
the results of equation 3.28 are obtained. 


vc (t) evsr(t-p) = vsr(t)vsr(t-p) (eqn 3.28) 
What ts:© Of interest is the average value of 


VST (t)vsr(t—p). The autocorrelation function (ACF) will be 
defined as that average value of vsr(t)vsr(t-p). The 
mathematical representation of the autocorrelation function 
for a continuous voltage is represented in equation 3.29 
while the autocorrelation function for a puise voltage is 


presented as equation 3.36. 


" 
ACF = Rvv(p) = 2 [ v(t) v(t-p) at (eqn 3.29) 
00 i 
ACF = Cvv(p) = [ v(t) v(t-p) at (eqn 3.30) 
— co 


The autocorrelation function is a measure of the degree 
to which two identical signals which are corrupted in some 
Manner are alike. The ACF of two signals. which are 


identical, not distorted in any way and occurring at exactly 


43 


the same time, 1.e., p = 0 has a maximum value. On the 
other hand there will be very little correlation between two 
identical time variant signals when the time difference 
between them is great. In this case the autocorrelation 
function is near zero. 

A Similar concept is used within the context of signal 
comparison. This concept is the crosscorrelation which is a 
measure of the degree to which 2 different signals are 
alike. When the crosscorrelation between voltage vsr and 
random noise voltage vn is 0 there is no relationship or 
correlation. 


Again let us assume that the signal at the receiver is a 


noise disrupted version of the Signal originally 
transmitted. 
vr(t) = vsr(t) + vn (t) (eqn 3.31) 


In the receiver , the time delayed version of the signal is 
applied to the incoming Signal and the average is formed as 


before and shown in Figure 3.17. 


Vr (+) = 
vsr(t)+r vn(+) 


| 

VSr(t)vsr (4 - 
—+ 
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Figure 3.17 Voltage Correlator 
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{[vsxr(t) + vn(t) Je{[ vsr(t-p) + vn(t-p) ] = (eqn 3.32) 
vsr (t)vsr(t-p) + vn(t)vsr(t-p) + 
vsr (t)vn(t-p) + vwn(t)vn(t-p) = 


Rvsr vsr(p) + Rvn vsr(p) + Rvsr vn(p) + Rvn vn (p) 


Rvn vsr(p) = 0 
Rvscr vn(p) = 0; because average vn = 0 
Therefore the result is simply Rvsr = vsr(p) + Rvn vn(p) 


or the autocorrelation function of the desired signal plus 
the autocorrelation function of the noise. Assuming the 
shape of the autocorrelation function of both components 
(both the original signal and the noise which vary with p, 
time)is known, what is done in practice is to vary the value 
of p until the ratio cf Rvsr vsr(p) / Rvn vn(p) iS a maxinun 
or the autocorrelation function of the signal is a maximum 
while the autocorrelation of the noise iS minimum and tne 


signal is recovered. 


De. THE MATCHED FILTER 


The matched filter is a linear filter which has the 
characteristic response desired for optimum reception of the 
desired signal. In other words, we desire the response of 
the system to the linear filter to be in some way 
proportional to the autocorrelation function of the desired 
Signal and in no way related to the autocorrelation function 
of the noise. This system could be illustrated as in Figure 
32,18 - : 
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Figure 3.18 Matched Filter Block Representation 


The question is what should the makeup or h(t), (1i.e., 
the response of the linear filter) be? The output of a 
. linear system is a convolution so the output of the linear 
filter, h(t) to an input vsr (t) would be vsr(t) * b fei 
where * denotes convolution. 

In the previous section it was stated that the 


autocorrelation function, Rvsr vsr(p), a function of p, is 


given by 
CO 
RVSI vsr(p) = [vsr (t)vsx(t-p) dt (eqn 3.33) 
—0o 


By simply changing variables, equation 3.33 can be rewritten 
as a £unction of @tine, +t: 
Co 
Rvsr vsr {t) = [ vst (p) sx (p-t) dp (eqn 3.34) 
-oCo 
From convoluticn theory it is known that the response of a 
linear filter as above, vsr(t) * h(t) can be expressed 
mathematically as in equation 3.35. 
GO 


vsr(t) * h(t) = [vsr(p) h(t-p) dp (eqn 3.35) 
—00 
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As stated earlier, the desired output of the linear 
filter is the autocorrelation function Rvsr vsr(t). Note 
the similarities between the last two equations. TE tan 
equation 3.35 a simple substitution of h(t) = vsr(-t) is 
performed, the desired results are obtained. 


(oe) 


vsr(t) * h(t) = [ vsr(p) vsr (p-t) dp (eqn 3.36) 
— 0d 
= Rvscr vsr(t) 
Thexcefore, for every signal that is desired to be 


recovered, the matched filter will do the job very nicely if 
the response of that filter, h(t) is equal to the inverse of 


the signal that is desired to be detected. 
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IV. ANALOG TO DIGITAL CONVERSION 

Key to the understanding of digital communications is 
how the analog information is converted to digital 
information. AS illustrated earlier, each analcg voltage or 
Signal can be represented by appropriately spaced sample 
values. These Sample values are still analog in that they 
can take on any value in the range of the original analog 
Signal. What is desired is to change this infinite set of 
decimal numbers toa finite set of decimal numbers. This is 
calied analog to digital conversion. Common types of A to D 
conversion include Pulse Code Modulation (PCM), Differential 
Pulse Code Modulation (DPCM), Delta Modulation (DM), Pulse 
Amplitude Modulation (PAM), Pulse Duration Modulation ~ 
(PDM),and Pulse Position Modulation (PPM). PCM, DPCM, and 
DM will ke examined in this chapter because of their 
widespread usage and easy application to digital technology. 


A. PULSE CODE MODULATION (PCM) 


PCM is the most widely used A to D conversion technigue. 
It involves assigning the sample value obtained in sampling 
to one- of several quantization levels within the tfrange of 
the voltage sampled and then representing those quantization 
levels by various binary code words. 

For example, the highest significant frequency component 
in human speech is 3300 Hz. In order to capture accurately 
the signal produced in speech, one would require a sampling 
rate greater than “or equdi “to” 238 or ~~ 2(3300) = 6600 
samples/sec. The telephone company uses 8000 samples/sec. 

B ="3300 “Hz 
f > 2B = 6600 samples/sec 
u 


se £ = 8000 samples/sec to ensure no aliasing 
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If a signal is sampled at the rate of 8000 samples/sec, then 
in the transmission of those samples, assuming the samples 
are transmitted immediately and not delayed, the 
transmission rate must be 

Ts = 1/f 

Ts 


Ts = 125 microsec/sample 


1/8000 samples/sec 


The number of quantization levels employed to represent 
the various analog sample values is arbitrary. The greater 
the number of levels, the more accurate the reconstruction 
of the originai signal but the more rapid the data 
transmission rate must be(in all cases there will be some 
error present after recovery). Consider again the example 


of voice as illustrated in Figure 4.1. 


v(r) 


Sample every 
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Figure 4.1 Samples of Voltage 


Assume it is desired to represent each sample value (decimal 
numoer) with an 8 bit binary code word. Then the number of 
quaatization levels is given by equation 4.1. 

The spacing between quantization levels is the range of 
voltages to be represented divided by the number of 


quaatization levels. The analog to digital conversion takes 
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Figure 4.2 Analog Nature of Samples 


#guantization levels = 2**n (eqn 4.1) 


for n = 8 bits/binary code word 
#quantization levels = 2**8 
= 256 


piace by determining into which quantization level the 
respective sample values fall and assigning that sample 
value the binary number associated with that quantization 
level. 

Now each sample value is assigned to an eight kit binary 
cod2 word. For voice transmission it has been shown that 
the minimum transmission rate was Ts = 125 microsec/sampie. 
The bit rate is then simply equal to the number of bits per 
sample times the inverse of the transmission rate as 
illustrated in equations 4.2 and 4.3. 

If voice was to be time multiplexed onto a single 
Channel, the bit rate for that channel would be the number 
of signals per channel times the bit rate for the most tiaie 


intensive signal as illustrated in equation 4.4. 
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bit rate = (n bits/sample) (1 sample/T sec) (eqn 4.2) 
bit rate = (8 bits/sample) (1 sample/125 microsec) 

= 64 kbps 

bit duration = 1/bit rate (eqn 4.3) 


= 1/64 kbps = 15.6 microsec/bit 


(#signals/channel) (bit rate/signal) = (eqn 4.4) 
bit rate/channel 


In other words, if 24 voice signals are to be 


multiplexed onto a single channel, the bit rate must be 


bit rate/channel 2 (eqn 4.5) 
(24 signals/channel) (64 kbps/signal) 
= 1536 Mbps 


As can be seen, the data transmission rate increases 
Significantly, necessitating better and better hardware and 
a wider and wider bandwidth. Advantages should be readily 
apparent for A to D conversion technigues or modulation 
techniques that reduce the bit rate required. 

Recovery of the analog signal is by a process called 
Digital to Analog Conversion. All that will be said about D 
to A conversion is that it is the inverse process of A to D 
conversion and that a slight error is always introduced 
during the process due to the discrete nature of the 


quantization process during A to D and D to A conversion. 


Be. DIFFERENTIAL PULSE CODE MODULATION (DPCM) 


In DPCM, what is converted to binary numbers is not the 
quantization level (decimal number) which each sample value 
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is represented by but the successive differences between 
quantization levels. The idea is that the range of maximum 
difference values will be smailer than the range of actual 
sample values. It is therefore possible to represent that 
range of delta values with fewer bits/binary word and 
therefore fewer quantization levels and ultimately a lower 
bit rate. [Ref. 11] | 

In actual practice, in a DPCM conversion technigue, it 
is a statistical estimate of each successive sampie vaiue 
which is subtracted from the actual value that is converted 
to binary code words. The result is the same in that the 
range of amplitudes is reduced and therefore fewer bits/word 
are reguired to represent the sample thereby lowering the 
data transmission rate required to transmit the signal. 
(Ref. 11} 


C. DELTA MODULATION (DM) 


Delta modulation is a form of DPCM where successive 
quantization levels of the output differ by only 1 bit. 
That is to say that any successive quantization level can be 
represented by varying only one bit of successive output 
binary code words. In other words a type of gray code is 
employed. This A to D technigue is implemented through the 
use of a DM coder or linear delta modulator. This DM coder 
approximates a given input Signal with a series of linear 
segments of uniform slope. A comparison is made between the 
value of this approximation andthe input signal at each 
sample increment. The sign of this difference value is what 
is encoded andis used to increment the DM coder in the 
direction of the input Signal. By using the differential 
sigh’ value of the input signal and the incremented 
approximation from the DM coder the linear approximation 
from the linear delta modulator is said to "track" the input 


Bz 


Signal. [Ref. 8] Slope overlcad of this type of modulation 
technique occurs when the slope of the incoming signal 
exceeds the ability of the DM system to follow the source at 


the sampling rate being utilized. [Ref. 11]. 
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Digital signal modulation techniques are the methods of 
encoding information for transmission utilizing digital 
technology. Factors affecting digital modulation include, 
but are not limited to, the physics of the method, hardware 
requirements, bandwidth considerations, power requirements, 
data transmission rates and error probabilities. ltr sa5 
these aspects which will be explored in the following 


secticns. 


A. DIGITAL HODULATION FORMATS 


Modulation is the technique by which the characteristics 
ot one waveform (called the carrier) are varied or modified 
by the characteristics of another (called the source). The 
Carrier waveform of interest is the sinusoid. It should be 
obvious that the attributes of a particular sinusoid that 
differentiate it from every other sinusoid are its 
amplitude, phase and frequency. It .follows..thateseae 
characteristics of the carrier waveform to be varied by the 
source are its amplitude, phase, frequency or any 
combination of the three. This gives rise to the broad 
general formats of Amplitude Shift Keying (ASK), Phase Shift 
Keying (PSK) and Freguency Shift Keying (FSK). All. other 
digital modulation technigues are variations of or 
combinations of these basic formats. Other factors involved 
in the description of the digital modulation technigue beiag 
employed inciude the number of bits being encoded at one 
time, the employment of error correction technigues and the 


baseband (source) waveform. 
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When each bit of the baseband waveform is individually 
encoded by any of the techniques previously mentioned (ASK, 
PSK, FSK), the technique being utilized is referred to as 
binary encoding. When more than 1 bit of the source code is 
modulated onto the carrier at one time it is called block 
encoding. Block encoding allows for one of nm = 2**k 
waveforms where k = number of bits. 

This paper wiil not go deeply into signal detection or 
demodulation techniques, however, a basic understanding of 
what is involved is necessary and again further defines the 
digital rodulation format being employed. In general terms, 
Signal detection is referred to as either coherent or 
noncoherent detection of the transmitted signal. Coherent 
detection, perhaps the easiest to understand, is when all 
possible waveforms of the modulated carrier waveform are 
available at the receiver and the waveform at the receiver 
is in phase with the transmitted carrier. Noncoherent 
detection is involved when the receiver does not have 
knowledge of the phase of the transmitted information and 
one of a number of phase estimation techniques must be 


employed for Signal recovery. 


B. HARDWARE 


Although significant advances have been made in recent 
years to improve the quality of hardware associated with 
satellite communications, most components are not what would 
be considered “otf ) the shevf items*. The hardware 
components most commonly referred to with regard to digital 
communications include the sampler, encoder, nodulator, 
Multiplexer and transmitter. There are variations of the 
above hardware requirements necessary for certain types of 
digital formatting, however, those exceptions will be 
addressed separately when the individual modulation 


technigues are examined. 
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Figure 5.1 Basic Hardware Component Block Diagram 


1. Sampler 


The sampler is the device or component which samples 
or extracts those single characteristic values ofeda 
haturally occurring analog signal which are ultimately 
encoded into a digital word by the encoder. The frequency 
at which this sampler must operate was derived in the 
discussion of the Sampling Theoren. It was determined that 
the Sampler must operate at a frequency greater than 2 times 
the highest significant Fourier frequency component of the 
source voltage. Modern solid state devices capable of 
taking thousands of samples/sec are available at modest 


costs. 


Z. “Encoder 


Often referred to as the A to D Converter (analog to 
digital converter), the encoder transforms the samples of 
the analog signal derived from the sampler into a digital 
format through one of the analog to digital conversion 
techniques described in the chapter on A to D conversion. 


These digital bits can be stored for later use, coded, 
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delayed or used immediately either individually or in groups 
to modify one of the characteristic qualities of the carrier 
waveforn. When sample values of the analog source signal 
are ccnverted into digital bits of information, they are 
Simply that, digital bits of information. The carrier can 
only be modified to represent the source information by 
interaction with another voitage. Therefore the value of 
the digital bit (binary), either 0 or 1, is used to generate 
a basetand waveform or voltage which does the actual 


modulation of the carrier waveforn. 
ae Common Baseband Waveforms 


It should be obvious that since it is desired to 
represent binary code words with a representative kaseband 
waveform what is required is two levels of voltage. There 
are two basic logic schemes for representation of the 
baseband waveforn. They are bipolar or unipolar logic. 
Bipolar logic involves representing the 0's or 1's of the 


birary codeword as either +V or -V as illustrated in Figure 


SEE aan 
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Figure 5.2 Bipolar Logic 


ea | 


The 0 or the 1 of the binary code word can be represented 
either as +V, -V or -V, +V respectively. This is a matter 
of convention but must be clearly understood in the various 
component designs in order to ensure compatibility between 
parts of the system. 

Unipolar logic utilizes a voltage to represent 
either of the possible binary digits and the absence of 
voltage to represent the other. Two common conventions of 


unipolar logic are represented in Figures 5.3 and 5.4. 
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Figure 5.3 Unipolar Positive Logic 


Again, which form of unipolar logic might be employed in the 
A to D converter is a matter of convention. 

Variations of these two basic encoding formats 
have the advantages of ease of generation and improved 
decoding and clock recoverability. Two variations commonly 
used in satellite communications are the Non Return to Zero 
(NRZ) and the Manchester waveforms. The NRZ wavetorms is 
Simply the voltage representation which corresfonds to the 
stream of bits represented in the logic scheme chosen. 
There is no transition as long as the same bit is present. 
Choosing bipolar logic as an example, a NRZ waveform can be 
described schematically as in Figure 5.5. 
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Figure 5.4 Unipolar Negative Logic 


v(+) vit)= eV > osts Ty 
Th= bit duration 
vit)=0 : elsewhere 


+V 


+ 


[ap np | 


in ee ee ee ee ee 


Figure 5.5 Non Return to Zero Waveforn 


The Manchester waveform of the same source 
voltage is represented with a transition at the midpoint of 
the bit duration from either +V¥V to -V or -vV to +#V. This 
form of coding has the advantage of clock synchronization in 
all cases of digital encoding. Schematically, the 


Manchester code can be represented as in Figure 5.6. 
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Figure 5.6 Manchester Waveforn 


3. Modulator 


The modulator in a digital communications system 1s 
nore commonly referred to as a noden 
(modulator/demodulator). This device does the actual 
transformation of the digital information into a waveform 
that can be transmitted from one point to another. # The 
modulation as mentioned earlier can modify either the 
amplitude, phase or frequency of the carrier wave. In 
Simple terms, the digital information derived from the 
analog source is mapped into the carrier wave. It is thas 
Mapping ‘that determines the power and bandwidth 
Characteristics of the moden. 


4. Muitiplexer 


Multiplexing is a common method for increasing the 
utility of a given communications channel. Variations in 
multiplexing technigues give rise to the multiple access 


techniques employed in satellite communications. 
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ae FPrequency-Division Multiplexing (FDM) 


Frequency-Division Multiplexing is a technique 
whereby the capacity of a communications channel is 
increased by adding one signal to another. These signais 
occupy discrete nonoverlapping frequency bands. A specific 
Signal is recovered by use of a band pass fiiter for the 
frequency band desired. [Ref. 7] 


b. Time-Division Multiplexing (TDM) 


Time-Division Multiplexing is a technigue used 
to increase the capacity of a given communication channel by 
adding signals in time. That. is, specific signals or 
portions of Signals are allocated specific time slots or 
portions of the carrier wave. These time allocations cannot 
overlap and the signal 1s recovered through proper 
synchronization of the recovery hardware with the 
multiplexer. [Ref. 7] 


c. Code-Division Multiplexing (CDM) 


Code-Division Multiplexing is a technique for 
increasing the capacity of a given communications channel 
through the assignment of a characteristic code to the 
digital signal before it is multiplexed in the time or 
frequency domain. Since the code used in multiplexing is 
Known at the demultiplexer, recovery of the digital source 
code could be accomplished through the use of a matched 


EEebeer: OF Correllator..i-[ Ref.) 7.) 


C. BANDWIDTH 


For purposes of discussion and for uniformity, when 
examining the various digital signal modulation techniques 
of interest, the bandwidth wili be defined as the highest 


Significant Fourier frequency component of the signal in 
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question. For digital communications, it can be shown that 
by the above definition, the bandwidth of the baseband 
waveform is 1/Tb, where Tb is the bit duration. For 
example: 
1. Voice sampled at the rate of 8000 samples/sec 
2. Each sample represented by an 8 bit ccde word 
3. Rb = bit rate = (8000 samples/sec) (8 bits/sample) 
Rb = 64 kbps 
4. Tb = bit duration = 1/Rb 
Tb = 15.6 microsec 
The Fcurier transform of the baseband waveform of duration 


15.6 microsec is represented in Figure 5.7. 
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Figure 5.7 Fourier Transform of Baseband Waveform 


It should be readily apparent that as the bit rate of the 
digital baseband waveform increases, so does the bandwidth. 
It is desireable to keep the bandwidth of the signal toa 
minimum due _ to Llianitatrons on availability of the 
electromagnetic spectrum. In addition, interference with 
Signals which occupy adjacent frequency bands is reduced and 
propogation limitations of the medium and hardware 
‘limitations such as self-inductance and capacitance are 


minimized. 
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D. SPECIFIC TECHNIQUES 


Now with some of the basics of digital signals already 
introduced it is possible to investigate some of the more 
significant modulation formats employed in digital satellite 
communications. The organization of this section will be 
around the three general types of modulation, i.e., phase, 
amplitude and frequency. Within each broad category, 
several techniques will be described under the subareas of 
binary encoding, where the information is modulated onto the 
Carrier bit by bit, and block encoding, where groups of bits 
do the modulation. Variations of these areas will be 
pointed out. The characteristics of the modulation 


technique will be introduced including, where practical, an 


analytical description, phasor representation, error 
peobability, spectral efficiency, advantages and 
disadvantages. 


1. Phase Shift Keying (PSK) 


Phase shift keying is a generic form of signal 
modulation where the baseband waveform is used to modify or 
Change the phase of the carrier. Phase, as stated earlier 
is one of the characteristics which distinguish one sinusoid 
from another. In general, phase shift keying offers the 


advantages of power and bandwidth efficiency. 
a. Binary Phase Shift Keying (BPSK) 


As the name suggests, Binary Phase Shift Keying 
(BPSK) results in a mcdulated carrier waveform consisting of 
two phases of the same sinusoid. The baseband waveform is 
used to change the phase of the carrier bit by bit from one 
phase to the other. In the case of bipolar logic, the bit 
SrEeCam Consists of +1. The general waveform of the carrier 


can be represented as a coSine wave of amplitude A, angular 
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frequency w and initial phase angle delta as in equation 
Seite 


vc(t) = A cos(wt + delta) (egn 5.1) 


Since two phases are represented in vc(t), it is 
customary to let them differ by pi radians. This can be 
represented by the mcdificaticn of equation 5.1 to account 
for a pi phase shift depending on whether the bit, v(t) to 


be aodulated is either a +1 or -1 as Shown in equation 5.2. 
vc (t) =A cos(wt + deita + pr/Zfi-yv (t)7) (egn 5.2) 


Equation 5.2 can be expanded to resuit in a simplified form 
Showing that the modulated carrier waveform is simply a 
product of the baseband waveform and the carrier. See 


equation 5.3 where v(t) represents the baseband waveforn. 
vco(t) = A v(t) cos (wt * detta) (egn 5.3) 


A phasor diagram can be constructed to represent 
the modulated carrier waveform in BPSK and is shown in 
Figure 5.8. AS can be seen by the phasor diagram in Figure 
5.8, the phase of the modulated waveform defends on the 
value of the baseband waveform and differs by pi depending 
on whether the value of the modulated bit is +1. 

Remembering the section on the properties of the 
Fourier transform, it iS possible to analyze the spectral 
and power efficiency of BPSK. Recall BPSK can be created by 
multiplication. Multiplication of two voltages is covered 
by the frequency translation property of the Fourier 
transform. Frequency translation results in double the 
bandwidth or spectral requirement of the translated waveform 


and half the power. In the case of a NRZ baseband waveforag, 
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Figure 5.8 Phasor Diagram of BPSK 


the bandwidth of the non-translated pulse, determined from 
the largest significant Fourier frequency component is 1/Tb, 
maere. Th is the bit duraticn. The bandwidth of the 
translated voltage (BPSK carrier) is 2/Tb. The power in a 
Sinuscid is A2/2 and none of this power capacity is lost in 
BPSK Since only the phase. cf. the original sinusoid is 
changed. For a Manchester baseband. waveforn, the 
non-tranSlated signal has bandwidth 2/Tb and the translated 
wave has bandwidth 4/Tb. 

The advantage of BPSK is that it is relatively 
efficient in bandwidth utilization and power capacity. Et 
does have the disadvantage of fairly large sidelobe 
components which contribute to interference with adjacent 
frequencies. This usualiy necessitates some type of 


filtering before the signal can be transmitted. 
b. Differential Binary Phase Shift Keying (DBPSK) 


The recovery of a Binary Phase Shift Keyed 
Signal requires the use of a phase coherent reference signal 
as described in the previous sections on autocorrelation and 


matched filters. This may not always be practical or 
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possitle soa technigue called Differential Binary Phase 
Shift Keying has be developed. DBPSK uses the same carrier 
type as BPSK and offers many of the same advantages and 
disadvantages. As stated, the primary difference is the 
lack of necessity for a phase coherent reference. 

The baseband waveform, v(t), for DBPSK as 
generated by comparing the phase of the next bit to be 
transmitted to that of the previously transmitted bit. For 
example, if a +1is encoded onto the carrier as v(t), it 
represents a particular phase cf the carrier. If the next 
bit to be sent is aiso a +1, there is no phase change, i.e., 
the bit remains the same over two successive bit durations. 
rf, however, the second bit to be sent isa -1, this 
represents a pi phase change and is decoded as a bit change 
over successive bit durations. This can be expressed 
mathematically as in equation 5.4 where bk represents the 
bit which wili be encoded on the carrier and bk-1 represents 
the previous bit already encoded and ak represents the 


present bit. 
bk = bEk-1 e ak (eqn 5.4) 


As can. be seen when bk = -1, there is a phase 
change and the bit changes from one bit to the next 
depending on the logic type being employed. If bk. =seae 
this represents the phase remaining constant over two 
successive bit duraticns, i.e., no bit change. See Figure 
5.9 [Ret 6. Recovery is effected by correlation of the 
received bit to the previously received bit. 

Differential Binary Phase Shift Keying (DBPSK) 
differs from Differentially Encoded Binary Phase Shift 
Keying (DEBPSK) in the recovery step ouly. DEPSK takes 
advantage of the modulation of phase shift in the recovery 


process while DEBPSK still utilizes a phase coherent carrier 
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Figure 5.9 Differential Binary Phase Shift Keying 
& . 
Signal in demodulation. Both cases use the same Scheme for 


differentiai encoding. 


c. Orthogonal Binary Phase Shift Keying (Orthogonal 
BPSK) 


In Binary Phase Shift Keying there is always a 
probability, with the interjection of noise on the 
transmitted signal that an individuai bit will be decoded 
incorrectly from that which was actually transmitted. In 
some environments, this probability is so high that special 
encoding technigues must be employed in order to reduce the 
probability that a bit or a binary code word will ke decoded 
hbo Crror. One such technigue is Orthogonal BPSK. 
Orthogonal BPSK is a type of group encoding scheme where 
groups or blocks of bits from the original signal are 
asSigned redundant bats? according to a predetermined 
assignment routine. This new expanded sequence of bits is 


then transmitted as before in BPSK. 
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For example, assume that it has been determined 
that the respective values of the original signal will be 
represented by a k bit binary code word. It has been 
demonstrated that there exists 2**k binary code words, k 
bits long in this type of arrangement. With Orthogonal BPSK 
there also exists 2**k orthogonal sequences of bits (called 
channel symbols, bauds or chips) [Ref. 6], which represent 


the k binary code words. However, the number of bits or 
chips in the sequence iS increased to n= 2**k, thereby 
providing the desired redundancy. Each of the “235 


sequences of hn chips 1S constructed to be orthogonal to 
every other sequence in order to assure maximum separation. 
See Figure 5.10 [Ref. 6]. 
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Figure 5.10 Orthogonal Binary Phase Shift Keying 


Decoding is aiso done by blocks through a bank 
or 2**k correlators. Since the decoding is dore in blocks, 
the prokability that the entire transmitted code word will 
be incorrectly decoded is reduced even though individual 
bits are incorrectly decoded. This is the major advantage 
of Orthogonal BPSK. Figure 5.11 [Ref. 6], shows the 
reduction in error probability with Orthogonal BPSK 
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Figure 5.11 Probability of Error for Orthogonal BPSK 


The major disadvantage to Orthogonal BPSK is 
enat. 2t results in a reduced data bit rate ata given 
transmission rate or the necessity for a higher transmission 
rate if the same bit rate is to be maintained. For example, 
suppose the bit rate of a PSK signal is R bits/sec. If each 
of the k bits in the PSK signal are represented by 2¥**k 


chips as is the case in Orthogonal BPSK modulation, then the 
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bit rate must be (2**k/k)R if the same data rate is to be 
maintained. Additionally, since the necessary bit rate 
increases, the bit duration decreases. It was shown in the 
secticn on bandwidth that as bit duration decreases as the 
result of a higher bit rate, an increase in bandwidth is the 
result. Therefore, the tradeoff in decreased transmission 
error comes at the expense of bandwidth and FEit rate. 
(Ref. 6] 


d. Quadrature Phase Shift Keying (QPSK) 


Quadrature Phase Shift Keying derives its name 
from a four-phase modulation of the carrier waveform. These 
modulations are achieved by simultaneously modulating the 
inputs from two bit streams onto a singie carrier. The two 
bit sequences could be from two separate sources or 
successive bits from a single source. To take advantage of 
the orthogonality of the sine and cosine functions, QPSK can 
be viewed as in eguation 5.5 as the sum of a BPSK modulated 
Sine and cosine. 


ve(t) = A vi(t) cos(wt + delta) + (eqn 5.5) 
A vz(t) sin(wt + delta) 


The signal can also be represented as in 
equation 5.6 where the phase angle, theta, of the modulated 


carrier is shown in eguation 5.7. 
vce(t) = 2**.5 A cos(wt + theta +t delta) (egn 5.6) 
theta = arctan (vi(t)/v2(t)) (eqn 5.7) 
The phases possible as a result of substitution of the 
possitle bits in birolar logic in equation 5.7 are? 28a 


degrees and +135 degrees. See Figure 5.12. 
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Figure 5.12 Modulation and Phases of QPSK 


Since two bits are being encoded at one time for 
transmission, the bit rate in QPSK is twice that of BPSK. 
The BPSK bit rate is equal to the QPSK symbol rate, (symbol 
rate being the time the guadrature bits exist on the 
Carrier). It is the symbol rate in QPSK which determines 
the bandwidth and since symbci rate ‘is equal to bit rate in 
BPSK, the bandwidth of the two modulation technigues is 
identical. The great advantage of QPSK lies in the 
bandwidth utilization. Two bits —of information are 
transmitted thereby effectively doubling the bit rate at no 
additional cost in bandwidth. A phasor trepresentation of 
QPSK is illustrated in Figure 5.13. Each quadrature 
component of the modulated signal contains half the power of 
the tctal carrier or (A2/2)/2 = A2/4. 


e. Offset Quadrature Phase Shift Keying (OQPSK) 


In Quadrature Phase Shift Keying, koth bits of 
the two baseband waveforms, which are to be modulated onto 
the carrier, change at the beginning of the respective 
symbol time. This allows the phase of the carrier to change 


up to 180 degrees for each symbol as iilustrated in the 
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Figure 5.13 Phasor Diagram of QPSK 


phasor diagram in Figure 5.13. The phase change of the 
modulated Signal can be limited to 90 degrees through a 
modulation technique called Offset Quadrature Phase Shift 
Keying. This iS accomplished by delaying the application of 
the second baseband channel for Ts/2 seconds, where Ts is 
Symbol duration as illustrated in Figure 5.14. Modulation 


is then accomplished as before in QPSK. 
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Figure 5.14 Offset Quadrature Phase Shift Keying 


W2 


The result is a QPSK modulated signal in which 
the maximum phase shift between successive bits is 90 
degrees. The bandwidth and power spectra of an OQPSK 
modulated signal are the same as that of a QPSK modulated 
signal. The advantage of OQESK is in the limit which is 
placed cn phase shift during encoding which simplifies the 
encoding hardware. Additionally OQPSK has spectral and 


interference advantages during decoding [Ref. 6]. 
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Figure 5.15 Phasor Diagram of OQPSK 


f. Multiple Phase Shift Keying (MPSK) 


Another form of digital modulation is knownas 
Multiple or M-ary Phase Shift Keying (MPSK). Once again, as 
in Orthogonal Binary Phase Shift Keying, k data bits are 
transmitted. The k bits represent the word length and there 
are m = 2**k different binary words, k bits long, for this 
type of modulation. In MPSK, the binary code word is used 
to vary the phase of the carrier. There are therefore a = 
2**k different phases in MPSK. These sequences of binary 
digits need not be orthogonal as in Orthogonal Binary Phase 


Le 


Shift Keying. The length of the transmitted symbol is only 
k bits long as compared to 2¥**k bits for Orthogonal BPSK. 
The phases therefore are not all orthogonal leading to 
greater difficulty in decoding and a higher probability that 
an individual code word will be incorrectly decoded, 
especiaily as m becomes large. 

As in other types of Phase Shift Keying, MPSK 
can be represented as in equation 5.8 where theta = phase = 
(2p1/8) i: 1=0, Vy0nicR— Ve 


vc(t) = A cos (wt + theta + delta) (egn 5.8) 


An advantage of MPSK is that as long as the 
symbol transmission rate does not increase, the carrier 
bandwidth does not increase even as the size of the binary 
code word increases. The disadvantage, as stated earlier, 
has tc de with error rates in decoding. Decoding requires a 
phase coherent reference of considerable stability 
especiaily as m increases and respective phases become 
closer together. Practical. limts for this type Soe 
modulation is m = 8 or a3 bit binary code word due to the 


complexity of the enccding and decoding hardware. 
2- Amplitude Shift Keying (ASK) 


Amplitude Shift Keying, as the name implies is 
modulation of the carrier through variation of its amplitude 
due tc variations in the baseband waveform. The information 
in the baseband waveform is thereby imparted to the carrier 
through this modulaticn. Since the anplitude of the carrier 
varies between successive symbols or binary code words, the 
power. also varies. Recovery is by comparison of the 
transmitted power during each symbol to the possible om = 


2**k power levels. 
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a. Multiple Amplitude Shift Keying (MASK) 


MASK is atype of group encoding where k bits 
are combined into a Single waveform for transmission and 
subseguent recovery. Since each symbol or binary code word 
contains k bits, there are m = 2**k different symbols and 
conseguently m different amplitudes of the carrier possible 
during the symbol duration time. An assignment scheme which 
maps the binary code words into the om different amplitudes 
would be to simply space them Ai = A/mi volts apart, where A 
is the maximum amplitude and mi represents the decimal 
equivalent of the binary code word from i = 1 to 2**k. 

The analytic form of a MASK waveform is 
represented in equaticn 5.9 where all the variables are the 
Same a presented earlier except for Ai which is equal to one 
of i different amplitudes depending on the bit sequence to 


be transmitted. 
vc(t) = Ai cos(wt +. delta) (eqn 5.9) 


MASK is not very popular in satellite 


communications since it depends on a carrier of very stable 


amplitude. This is not generally practical under actual 
coniitions. However, MASK does have the spectrai advantage 
of a constant bandwidth of 2/Ts, Ts equal to symbol 


duration, even as the binary code word increases in length. 


1 


De QuUaGrature Amplitude Shift Keying (QASK) 


Quadrature Amplitude Shift Keying (QASK) isa 
hybrid digital signal modulaticn technique. It represents 
separate amplitude modulation of the quadrature components 
of a ccmmon carrier, In that sense it is both ASK and PSK. 
Implementation of QASK involves simultaneous application o£ 
M-ary Amplitude Shift Keying to the quadratures. See 


oS 


equation 5.10, where Al and A2 are derived according to an 


assignment scheme aS in MASK. 


vc(t) = Al cos(wt + delta) + (eqn 5.10) 
A2 sin(wt + delta) 


By this’- method, the effective bit rate is 
doubled over that of ordinary MASK. In other words, QASK 
results in the same data rate as modulating 2k bits onto the 
carrier at the same time or during the same symbol duration 
with MASK. The advantage to QASK lies in the fact that this 
increase in bit rate comes at no further increase in 
bandwidth. The disadvantages associated with MASK are still 
present with an additional increase in the complexity of the 
decoding equipment. Decoding of QASK must be done in two 
steps. First the signal is recovered in phase through a 
phase-coherent correlator and then each amplitude modulated 
Signal is recovered as before in MASK by comparison of.the 


power levels of the received signal. 


3. Frequency Shift Keying (FSK) 


Frequency Shift Keying (FSK), is the generic tera 
used to describe a number of digital signal modulation 
techniques which cause variations in the carrier frequency 
by interaction with the baseband waveforn. Decoding is 
generally through measurement of the freguency of the 


received signal. 


a. Mimimum Shirt Keying or fFast Frequency Shift 
Keying 


Power requirement for the transmission of 
retransmission of a satellite signal is extremely important. 


It is desireable to limit power required to accurately 
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transmit a piece of information to the minimum possible 
consistent with allowable error rates. For this reason, 
detection of FSK signals is limited to coherent methods as 
power required increases significantly for non-coherent 
methods. 

MSK or FFSK are identical and represent a 
modulation technique known as continuous phase FSK. They 
get their names from the fact that "fast" indicates more 
bits per second can be transmitted in a given bandwidth than 
ordinary BPSK and “"mininum" refers to the minimum modulation 
index for which orthogonal signalling occurs [Ref. 7}. 

Analytically, FFSK can be expressed as in 
equation 5.11 and equation 5.12. 


vc(t) = A cos((we + deltaw)t) (eqn 5.11) 


vce(t) = A cos(+deitaw t) cos(wc t) - (eqn 5.12) 


A sin(+deltaw t) sin(we t) 


AS is noted, FFSK can be envisioned as the 
summation cf separately modulated in-phase and quadrature 
components of the carrier by the baseband waveforn. The 
technique is Similar to that used in OQPSK. The frequency 
variation in the carrier waveform is between wc + deltaw and 
woe - deltaw. Maximum separation in phase of these two 
frequency components occurs at pi as noted in equation 5.13, 


where Tb represents the bit or symbol duration. 
(w2-w1)Tb = pi (eqn 5.13) 
This can be converted to the modulation index 


menticned above simply by converting to frequency as in 


equation 5.14. 


Rid) 


h -=(£2-£7 bee 425 (eqn 5.14) 


For FFSK, the frequency deviation or separation 
between frequencies of the transmitted carrier is exactly 
1/2Tb and the deviation from the carrier frequency is 1/4Tb. 
FFSK can therefore be rewritten from equation 5.12 as in 
equation 5.15 by substitution for deltaw = 2pi(1/4Tb). 


vc(t) = A cos(£pi £7275) cos(Zpi fe t) = (eqn 5.15) 
A sin(tpi t/2Tb) sin(2pi fc t) 


The error rate in this modulation technique is 
identical to that of BPSK. Frequency Shift Keying can be 
implemented to incorporate any given modulation index by 
simple calculation of the frequency separation with equation 
5.14 and: substitution |anto ” equation . Sats. For > \Amamy 
Frequency Shift Keying, the separation between frequencies 
must te at least 1/TsS, where Ts is the symbol duration, in 
order to avoid carrier energy from one frequency being 
incorrectly interpreted as carrier energy from another 
frequency. This results in a modulation index of i. Figure 
5.16 represents the phasor diagram of an FFSK modulated 


signal. 
4. Quadrature Partial Response Signalling (QPRS) 


Quadrature Fartial Response Signalling (QPRS) 
represents a specific type of the general class of digital 
signal modulation techniques called Partial Response 
Signalling. Partial Response Signalling is a method in 
which a controlled amount of intersymbol interference (ISI) 
is allowed during the encoding process. Previously any 
overlap in successive Signals resulted in aliasing and an 


increased probability of error in decoding. The idea was 
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Figure 5. 16 Phasor Diagram of FFSK or MSK 


first introduced by Lender in 1963 [Ref. 12]. His concept 
was that in knowing and controlling the amount of 
interference allowed during the encoding process, 
compensations can then be made for the ISI at the receiver. 
Allowing for a limited amount of ISI makes it possibie to 
transmit at rates eyual to the Nyquist Cate, something that 
is not possible with many other systems [Ref. 13]. 

QPRS iS implemented, as with previous quadrature 


systems, through simultaneous modulation of the quadrature 


componeats of a common carrier. However, this time the 
modulation is accomplished with a PRS systen. The 
modulation of the respective quadrature components 


represents the impulse response of one of a number of linear 
filters to the bits cf the baseband waveforn. The linear 
filter employed depends on the class of the Partiai Response 
Signalling system being used. There are 5 classes of linear 
filters commonly used in PRS and they are illustrated 
analytically in Table 1 and graphically in Figure 5.17. 
Output of a QPRS system can be represented as in 
equation 5.16, where an and bn represent the nth bit to be 


modulated and h(t - nTIs) represents the contribution of the 
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TABLE 1 
QOPRS SYSTEM POLYNOMIALS 


SYSTEM POLYNOMIAL} FREQUEN DY RKESPUASE IMPLLSe RESPONSE hit) 
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Figure 5.17 QPRS Linear Filter Impulse Responses 


impulse response of the linear filter at time t for bit n. 
[Ref. 14] 
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(ore) 
vg{t) = yan h(t - n Ts) cos(wt + delta) (eqn 5.16) 


NF-OO : 
) kn h(t - n Ts) sin(wt + delta) 
n>=0O 


QPRS has the advantage of rapid transmiSsion rates 
with no increase in bandwidth and excellent error handling 
performance. The cost of such performance improvement comes 
in the form of a higher required signal to noise (SNR) ratio 
when compared to other binary systems. 
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Due to the nature of digital signal modulation, itis 
possitle to simulate systems through the use of the digital 
computer. This chapter is a description of the computer 
Simulation included in the Appendix of this thesis. 


A. GENERAL DESCRIPTION OF THE PROGRAM 


The program was constructed using top-down design and 
modular programming. FORTRAN was selected as the 
programming language due to its capabilities in the area of 
humerical operations andthe availabiiity of mathematical 
routines to be used as modules in the progran. Testing was 
accomplished on each module as the program was developed. 
The program was designed to be used interactively in order 
to allow for instructional use. Although the program was 
made user friendly in that error checking is accomplished on 
user inputs, care must be taken to insure instructions are 
followed correctly. 

The program consists of twenty modules, sixteen of which . 
were written by the author, three which were taken from the 
Double precision International Mathematics and Statistics 
Library (IMSLDP) and one which was taken from a NON-IMSL 
library which resides on the IBM 3033 iocated at the Naval 
Postgraduate School. The IBM 3033 was used exclusively for 
development and testing. 

The development was accomplished using the WATFIV 
compiler; however, the program was written to operate with 
the VS FORTRAN compiler as well. Testing and operation was 
done using the VS FORTRAN compiler. Figure 6.1 is a bilock 
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diagram of the relationshifr of the modules for a 


representative modulation subroutine. 


B. MAIN CONTROL MODULE 


MAIN operates as the control module for the entire 
computer Simulation and accomplishes limited output. ze 
introduces the program and allows the user to input various 
parameters of the digital signal modulation technique to be 
Simulated and various other control functions. These inputs 
include: 

1. Modulation technique 

2. Class of QPRS system (when appropriate) 

3. Digital logic scheme 

4. Baud or symbol rate 

5. Bits/binary code word (when appropriate) 

6. Carrier frequency 

7. Number of samples to be generated 

8. Carrier max amplitude 

9. Initial phase angle 

10. Use or random Nhumber generator or no 

11. Seed for RNG (when appropriate) 

12. Number of repetitions of the simulation 

Once the uSer has input the desired characteristics of 
the modulation to be simulated, MAIN calls the appropriate 
module to begin the actual ‘calculation. The values 
necessary to calculate the time series of the signal are 
passed to the subroutines aS parameters. HAIN “-calis' the 
required subroutine the number of times the user specifies 
as repetitions. Each call to a subroutine produces’ the 
Discrete Fourier Transform and amplitude spectrum of the 
Signal. It is these respective amplitude spectrums which 
produce the statistics in the final output after MAIN calls 


the statistics subroutine the final time. 
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Figure 6.1 Representative Block Diagram 
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C. SUBROUTINE BPSK 


This description of SUBROUTINE BPSK, a module that 
Simulates Binary Phase Shift Keying, will serve as a guide 
on the construction and operation of subsequent modules 
involved in the actual signal modulation. A complete 
description including its interaction with other progran 
subroutines will be included. Since other signal modulation 
subroutines interact in the same manner within the frogran, 
subseguent descriptions will concentrate on differences to 
this basic module. 

SUBROUTINE BPSK begins with variable declarations as 
with all FORTRAN programs. All variables are declared to be 
either type integer or double precision. All calculations 
in the program are carried out in double precision accuracy. 

Variable initialization is next. Significant variable 
initialized at this point are TIME, the time of the first 
sample which is set to 0. STEP, the delta time between 
Samples is set to the normalized Nyquist sampling rate. 
OMEGA, the angular frequency and DELTA, the initial phase 
angle are computed in radians. NBAUD, a variable wiich 


keeps track of the number of symbols which have been 


nodulated, 1S ett toed worl, theca first symboloj as being 
modulated. #$NBAUD remains at its present value until total 
elapsed time exceed 1 symbol duration or BAUDD, the 


normalized baud rate. 

Subsequently, ani anatial)’:callbas «madeto:; SUBROUTINE 
STREAM to receive. the value of the first bit to be 
modulated, The modulation is then carried out according to 
equation 6.1. 

The values of each Sample are assigned to an array for 
storage along with the corresponding time increment. 
Modulation continues with the first drawn bit until 1 bit or 
symbol duration is exceeded when STREAM is calied to draw 
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vc(t) = A vit) cos(wt + delta) (eqn 6.1) 


A = anplitude 

v(t) = bit to be modulated 
w = angular frequency 

t = time 


delta = initial phase angle 


another bit. This entire process iS repeated until the 
number of samples specified by the user in MAIN is reached. 

At thas peiat, only on the first repetition of the 
Simulation, SUBROUTINE PLOT is called which gives the user 
the opportunity to plot the time series of the simulated 
Signal. Upon return from PLOT the Subroutine calls the IMSL 
routine FFTRC to generate the Discrete Fourier Traasform of 
the time series produced. This is accomplished on each 
repetition of the program. 

Again on the first repetition only, information on the 
nhumber.of the principle harmonic of the FFT is displayed. 
Plot is called again to give the user the opportunity to 
plot the amplitude spectrun. The values of the FFT and the 
amplitude spectrum are computed on each successive 
repetition of the program and the amplitude spectrum is 
added to the statistics being accumulated by a call to 
SUBROUTENE STATIS. 


D. SUBROUTINE DBPSK 


SUBROUTINE DBPSK simulated Differential Binary Phase 
Shift Keying and is essentially the same as SUBROUTINE BPSK. 
The difference lies in the fact that what is encoded during 
the modulation process is the product of the random bit 
drawn and the value of the bit previously modulated. A 
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reference bit equal to +1 is used to determine the first bit 
to be modulated. In this manner, modulation of +1 means no 
change occurs between successive bits and the phase of the 
Carrier remains the same. Modulation of a -1 indicates a 
change of bits has occurred and is indicated by a change of 
phase of the carrier between successive bit durations. The 
rest of the module remains unchanged. 


E. SUBROUTINE OBPSK 


Orthogonal Binary Phase Shift Keying is accomplished in 
this module. Again SUBROUTINE OBPSK is constructed 
essentially the same as SUBROUTINE BPSK. In Orthogonal BPSK 
the only bit streams that are modulated are n bits long 
representing each binary code word, where n = 2¥**k (k = 
humber of bits per word is limited to 6 bits in MAIN). Each 
sequence of n bits 1S orthogonal to every other allowed 
sequence. 

The way this is accomplished in this module is first to 
draw a series of k random bits from SUBROUTINE STREAM. This 
series of k 1's and 0's is then changed to its decimal 
eguivaient from 0 to 2**k-1 and saved for later use. 

The program continues with a call to SUBROUTINE ORTHO 
which generates an n X n orthogonal matrix of +1"s and -1's 
through the use of the Hadamard matrix and the Kronecker 
product [Ref. 15]. Tie! Hadamard *natrix* Ps’ a“ “2--K° 2 
orthogonal matrix of 1's and -1's. The Kronecker product is 
the matrix which is formed when a matrix is expanded to 
twice its original size by replacing each element of the 
Hadamard matrix by the product of the Hadamard matrix and 
the original matrix. The value of the decimal equivalent to 
the k random bits is used to identiry the row of the matrix 
and the respective columns represent the value of the binary 


digit which is modulated onto the carrier by BPSK. THis 
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sequence continues until either one bit duration is 
exceeded, when another orthogonal bit is received from ORTHO 
or when all the orthogonal bits in a row are modulated. Then 
another seguence of k random bits is converted to its 
decimal equivalent and a new row and column of the 
orthogonal matrix is identified. Execution continues in 
this manner until all required samples have been produced. 
It shculd be noted that the bit rate for Orthogonal BPSK is 
nh times the baud rate. This significantly increases the 
Signal bandwidth requirement and decreases the time 
increment between successive Sampies in the simulation. fThe 
rest of the module is the same as those previously 


presented. 


F. SUBRCUTINE QPSK 


Quadrature Phase Shift Keying, Simulated in SUBROUTINE 
QPSK, is accomplished through simultaneous BPSK modulation 
of the guadrature components of the carrier. Two randon 
bits are drawn from successive calls of SUBROUTINE STREAM 
and- each. bit is 2n turn modulated onto the carrier 
components and the sum formed. This process is repeated at 
time intervals equai to the normalized Nyquist sampling rate 
until TIME exceeds one symbol duration, BAUDD. At this time 
two new random bits are drawn and the modulation continues 
until all required samples are produced. The rest of the 
modules is the same as those before. Eguation 6.2 


represents the analytical expression used to sinulate 


G. SUBROUTINE OQPSK 


SUBROUTINE OQPSK is essentially the same as SUBROUTINE 


QPSK with these minor deviations. The time cf the first 
sample is artificially initialized as TIME = .5(BAUDD) of 
1/2 the symbol duration. In this manner, the two random 
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vc(t) = A v1(t) cos(wt + delta) + (eqn 6.2) 
A v2(t) sin (wt + delta) 


v1(t) = in-phase bit 


v2 (t) quadrature bit 


bits first being modulated are both known to have existed on 
the carrier for 1/2 the symbol duration. The first of thest 
two bits only remains on the carrier until TIME = 1(BAUDD) 


when a third random bit is drawn and modulation begins with 


Enaes bit. The second random bit is allowed to remain until 
TIME = 1.5(BAUDD) when a fourth random bit is drawn to 
Fepiace it and SO on. In this manner the two bits being 


modulated are offset in the time they change by 1/2 the 
symbol duration. This necessitates keeping track of the 
humber of bits modulated on both the in-phase and quadrature 
component of the carrier. The remaining portion of the 
module is the same as SUBROUTINE QPSK. 


H. SUBROUTINE MPSK 


This module simulates M-ary Phase Shift Keying. This is 
accomplished by changing the phase of the carrier to any one 
of n =2**k phases where k is the number of bits in the 
binary code word. The max length of the binary code word is 
limited to 10 inthe MAIN progran. This modulation is 
accomplished by use of equation 6.3 to simulate the MPSK 
systen. 

The program first draws k random bits and converts then 
to the decimal equivalent. At this point modulation begins 
at TIME = 0 according to equation 6.3. Modulation continues 
until one symbol duration is exceeded when a new set of k 


random bits is drawn, conversion to decimal takes place and 
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vc(t) = A cos{[wt + (2pi m)/n + delta] (eqn 6. 3) 


decimal equivalent of binary code word 


i) 


2**k; k = bits in binary code word 


modulation continues or until all desired samples have been 


produced. The program continues aS in previous subroutines. 


I. SUBROUTINE MASK 


The principles cf operation of SUBROUTINE MASK is 
Similar to SUBROUTINE MPSK which also involves a group 
encoding system. SUBROUTINE MASK allows for M-ary Amplitude 
Shift Keying. The analytic expression used in the 
Simulaticn iS provided in equation 6.4 Ai is one of 2¥**k 


egually spaced amplitudes. 
vc(t) = Ai cos(wt + delta) _ (egn 6.4) 


Once again a set of k random bits is drawn from 
SUBROUTINE STREAM. The decimal equivalent of the k bits is 
computed and used to adjust the amplitude Ai according to 
the asSignment routine expressed in equation 6.5. 


Ai 2 A/n (eqn 6.5) 
A = max anplitude 
w =.1 toby “nh = .2**k 
This process is repeated at each symbol duration until 
all desired samples have been produced. The rest of the 


module remains unchanged. 
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J. SUBROUTINE QASK 


SUBROUTINE QASK combines the elements of the previously 
described amplitude modulation technique with the now 
familiar quadrature modulation technique. Two separate 
streams of random numbers are generated and converted to 
their decimal fora. These numbers are used to produce two 
separate amplitudes with the same assignment scheme used in 
MASK. These amplitudes modulate the quadrature components 
of a common carrier and the output is formed by the sun of 
the quadratures. Modulation continues in the above manner 
until a symbol duration elapses when two new amplitudes are 
calculated. Processing terminates when all required samples 


are ccmputed. 


Ke. SUBROUTINE ASK 


Minimum Shift Keying is frequency shift keying in which 
the modulation index is 1/2. The modulation index is 


represented in equaticn 6.6 


h = .5 = (deltaf) Tb (eqn 6.6) 


deltaf = frequency separation 


Tb = bit duraticn 


The frequency deviation frem the carrier is .5(deltaf) = 
1/4Tb. Converting this frequency deviation to radians 
yields pi/2Tb. This angular frequency deviation from the 
central carrier is either +pi/1Tb depending on the value of 
the bit to be modulated. 

The simulation in SUBROUTINE MSK is accomplished by 
drawing a random bit and generating the signal according to 


equation 6.7. 
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vc (t) = A cos[{ (w +pi/2Tb)t + delta] {egn 6.7) 


Once again modulation continues until one bit duration 
is elapsed when another bit is drawn. Modulation stops when 
all reguired samples are produced. The rest of the 


subroutine remains unchanged. 


L. SUBROUTINE MFSK 


This subroutine modulates a signal simulating M-ary 


Frequency Shift Keying. The separation between adjacent 
freyuencies is estabiished as 1/Ts, where Ts is the symbol 
duration. This amounts to a moduiation index egual to 1. 


Initially the ertire range of frequency deviation from the 
Carrier is calculated as (n-1)Ts. The mean frequency 
deviation is then calcuiated. In other words the range of 
frequencies of the modulated signal is the frequency of the 
central carrier +1/2 the magnitude of the entire range of 
frequency deviation. The minimum frequency is used as the 
base for computing the frequency to be used to modulate the 
Signal during each respective symbol duration. 

At this point in the subroutine a series of k random 
bits (k<10) iS drawn and converted to decimal. This decinal 
number is multiplied by the frequency separation and added 
to the minimum frequency and converted to radians. The 
analyticai expressicn of the signal produced with this 


Subroutine is illustrated in eguation 6.8. 


ve(t) = Aicosf>€w to2pi emt) t?4? dettaq (eqn 6.8) 


B 
Fh 
iT} 


modulation frequency 


This modulation continues until a symbol duration 


elapses when a new frequency deviation or af is calculated 


a2 


or mcdulation stops do to completion of all required 


samples. 


H. SUBROUTINE QPRS 


This mcdule simulates 5 different classes of Quadrature 
Partial Response Signalling systems. The technique employed 
involves the summation, at ¢€ach time increment, of all 
responses to the linear fiiter representing the class. This 
is accomplished for all of the bits being modulated during 
the duration of the simulation. 

Tnitalily, two arrays of random bits are generated 
representing those bits which would be modulated onto the 
quadrature components of the carrier during the length of 
time the signal is to be simulated. The the impulse 
response to each of these n bits is calculated for the time 
of the respective sample. These impulse responses represent 
the modulation for the bit in guestion onto the respective 
porticns of the carrier and the sum is formed according to 


eguation 6.9. 


CO 
we (t) = >A h’'(® —*'nis) cosiwt + delta) + (eqn 6.9) 
CO N=~6O 
) A h(t - ats) sin (wt + delta) 


n= 


h(t - nIs) = impulse response of nth bit at time t 


The time of the first sample is initialized at 6(BAUDD) 
in order to ensure that when the first sample is taken, 
contributions from bits modulated at time 0 are also summed 
do to the overlapping nature of the QPR Signals. The value 
of the last sample is also formed by the sum of the 
guadrature components existing at 6({BAUDD) atrer it, is 
produced. The rest of the program operation fremains the 


Same aS previous modules. 
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NHN. SUBROUTINE STREAS 


This subroutine interacts with the IMSL random number 
generator SUBROUTINE GGUBFS to produce a random nunber 
between 0 and 1 and make assignment on the basis of the 
value of this random number to random bits according to the 
digital logic scheme specified in the MAIN program. It also 
enables the user to manually insert binary digits if that 


option was specified previously in MAIN. 


O. SUBROUTINE ORTHO 


CRETHO produces an nh Xn orthogonal matrix from the 
Hadamard matrix by forming successive Kronecker products 
{Ref. 15]. It also selects the appropriate row and column 
of the orthogonal bit to be soduiated and passes this bit 
back to SUBROUTINE OEFPSK. Foint of entry to SUBROUTINE 
ORTHO is controlled by a flag set in OBPSK. 


P. SUBROUTINE PLOT 


PIOT is an interactive module that allows the user to 
determine whether or not a graph of the output of program 
calculations is to be produced. The user is also able to 
selectively vary certain parameters associated with the 
plot. PLOT calls SUBROUTINE UTPLOT which actually produces 
the graph and performs the output. UTPLOT is a NON-IMSL 


library routine from the Naval Postgraduate School. 


Q- SUBROUTINE STATS 


The statistical calculations associated with the 
ampiitude spectrums of the various functions generated on 
successive repetitions of the program are computed in 
SUBROUTINE STATS. Upon completion of each repetition, each 


module calls STATS where the sum, sum of the squares, sum of 
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the cubes and sum of the gquartics of each element of the 
amplitude spectrum of the Signal produced are computed and 
stored. when the last repetition of the program is 
complete, MAIN instructs STATS to compute and output the 
fumal statistics;: i.e., mean, Variance, skewness and 
kurtosis of each component of tke amplitude spectrun 
according to the estimations contained in [Ref. 16]. 
Additionally the gean variance and variance of the variance 
is calculated and output. Point of entrance to this 
subroutine is controlled by flags set inthe individual 


modulation subroutines or in MAIN. 
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VII. STATISTICS OF THE FAST FOURIER TRANSFORMS 


A. DESCRIPTION OF THE PROBLEM 


As stated in Chapter II, if it can be shown that the 
statistics of the FFT can be somehow linked to a particular 
digital signal modulation technigue, then the modulation 
technigue can be identified on the basis of those statistics 
alone. Although it is beyond the scope of this thesis to 
actually derive those relationships if they exist, it seen 
prudent to attempt to determine if there is statistical 
differences between the components of the FFT's as they are 
derived in the enclosed computer simulations. 

The statistic chosen to do the hypothesis testing is the 
F-test since the true mean and variance of the distribution 
need not be known [Reft. 17]. In addition, the information 
necessary to calculate the F-statistic is readily 
accumulated in SUBROUTINE STATS. Calculation of the 
F-statistic necessitated the writing of a small computer 
program to be used for that purpose once the output from the 


Main frogram was generated and reformatted. 


B. ANALYSIS-OF-VARI ANCE 


If there were no differences attributable to the 
modulation technique, then it would be reasonable to assume 
that for a certain set of characteristics, the components of 
successive FFT's would be the same. In other words, if a 
Signal was modulated by BPSK and a statistically significant 
humber of FFT's were generated of the signal, then the 
statistics of those FFT's would be assumed to be related. 
On the other hand, if it can be shown that the statistics 


are not related (i.e€., not from the same distribution), then 
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significance can be placed in the variation among modulation 


technigues. 


1. The F-Distribution 


The assumptions necessary to use the F-distribution 

are: 

1. Normal distribution 

2. Random samples 

3. Independent samples 
These asSumptions are not too difficult to intuitively 
accept in tne model that will be proposed. Once again it 
would be expected that the FFI's of successive identically 
modulated signals to. be related. For the simulations 
conducted as part of this test, they were all generated 
using bits from a random number generator for which it can 
be shown that each seguence of bits passes statistical tests 
for randomness and independence. The assumption of 
normality could also be argued on the basis of the central 
limit theoren. 

The value of the F-statistic with which comparison 
will be made is 1.51 for 14 degrees of freedom in the 
humerator, an infinite number of degrees of freedom in the 
denominator and a 90% confidence region. How these values 
‘were oktained will be detailed under the design of the 


experiment. 
2. Design of the Experiment 


The experimental data used in the computation of the 
statistics was chosen to minimize the random contributions 
of variables other than the. modulation technigue. The 
computer simulations included in this thesis were used to 
generate the statistics and a separate program also included 
in the appendix was used to calculate the F-statistic. The 


variables in the experiment include: 
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1. Modulation technique 

2. Logic type 

3. Baad rate 

4. Bits per binary code word 

5. Carrier frequency 

6. Carrier amplitude 

7. Initial phase angle 

8. Number of samples generated 

9. Time between samples 

10. Seed for random number generator 

11. Number of repetitions or trials 

Fifteen different modulation techniques were 
compared. In all cases bipolar logic was simulated, the 
baud rate was held constant at 1200 baud and the maxinum 
carrier amplitude was established at 1 volt... Additionaily 
the phase angle of ail simulations was 0 degrees and the 
seed for the random number generator was 1 thereby ensuring 
the same random seguence of bits. The number of bits per 
binary code word did vary among the modulation techniques. 
When the modulation technique was M-ary, the bits per code 
word was always 3 so it is possible to infer that bits per 
code word 1s a function of the modulation technique. The 
time between sanples was the normalized Nyguist sampling 
rate. This did vary between modulation technigues but was 
necessary to derive an accurate FFT. Carrier frequency also 
varied tetween Simulations but was always twice the lowest 
carrier frequency recommended in the computer simulation. 
In most cases this was 2400 Hz, or twice the baud rate. 
Finally each simulation was repeated 100 times and the 
statistics of the FFT's based on those 100 repetitions. The 
number of samples produced was always 64 so there are 33 
components of the FFT. 
What this amounts to can be illustrated through an 

Analysis-of-Variance table as shown in Table 2 [Ref. 17]. 
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TABLE 2 
ANALYSIS-OF-VARIANCE TABLE 


n Observations in Each of r Groups 


| | Observations | Sum _| 


a cn a a cs a ca A a A A ces a es RS SO De ce ee ee ed 


reo yee rece | 


In the case of the Simulation described above, this amounts 
to 100 observations (repetitions or triais) from each of 15 
groups (modulation technigues) for each of the 33 components 
et the FFT. The degrees of freedom in the numerator is 
eqaal ta (15-1) or 14 while the degrees of freedcm in the 
denominator is equal to 15(100-1) or 1485. Table values for 
the F-distribution use infinity as the degree of freedom tor 
values greater than 120. 

An F-test was also performed to determine if there 
is a relationship among the mean variance, mean Skewness or 
mean kurtosis of the 15 different modulation techniques. 
Again the assumption is made that these statistics would be 
normaliy distributed from modulation technigues which had 
the same or statistically similar components of their FFT's. 
In this case there were again 15 modulation techniques to 
compare (14 degrees of freedom) but only 33 observations 
(the variance, skewness or kurtosis of the respective 
elements of the FFT). This number ot observations yields 


15(33-1) = 480 degrees of freedom in the denominator. 


og 


3. Hypothesis and Hypothesis Testing 


The hypothesis posed for the model is that the means 
of the individual ccmponents of the FFT's from the 15 
modulation technigues are fron the same distribution. 
Therefore, the respective means must be equal. Also tested 
is that the mean variance, mean skewness and mean kurtosis 
of all components of the FFT's of each modulation technique 
are equal. These hypothesis are illustrated in equations 
7.1. through © equation 7.4. Each hypothesis is tested and 
compared to the F-distribution selecting a rejection region 
or level —of--signpficance..o£).3s This equates to an 
F-statistic of 1.51 for 14 degrees of freedom in the 
Mumerator and an infinite number of degrees of freedom in 


the dencminator. 
4. Results 


The means of the 33 ccmponents of the FFT for each 
of the 15 modulation techniques were compared using the 
F-distribution. The results are shown in Table 3. Table 4 
shows the F-statistics associated with the comparison of the 
mean of the variance, mean of the skewness and mean of the 
kurtosis for the respective grours. In addition a complete 
summation of the results of the F-test are included in the 


appendix. 


C. CCNCLUSION 


The results of . the F-test indicate that the aull 
hypothesis must be rejected at the .9 significance level and 
the aiternate aypothesis accepted that the means are not 
egual for FFT components 9-14. In addition the F-statistic 
for the comparison of mean variance, mean Skewness and mean 


kurtosis all fall in the rejection area. 
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Ho: XBAR11 = XBAR21 = e @ e = XBARij (eqn 7.1) 


H1: XBAR11 # XBAR21# ee e # XBARij 


i = ith modulation technique 


= jth component of the FIT 


Hos MVAR1 = MVAR2 MVAR15 (eqn 7.2) 


Hi: MVAR1 # MVAR2 # © e @ # MVAR 15 


HO: MSKEW1 = MSKEW2 = e e = MSKEW15 (eqn 7.3) 
H1: MSKEW1 # MSKEW2 # e © © # MSKEW15 


Ho: MKUR1 = MKOUR2Z =e @ e 


MKUR15 (eqn 7.4) 


H1: MKUR1 # MKUR2 # e @© @ # MKUKI5 


This indicates that there is a difference between the 
Statistics of of the FFT's of the respective modulation 
technigue. Since all the parameters used in the generation 
of these statistics were held essentially constant among the 
trials, it can be assumed that these differences are due to 
the modulation technique emrfloyed. It has not been 
established yet what those differences may be. This is an 
area where future research and study is warranted. 

The results of this test would seem to point to those 
components of the FFT which offer the best opportunity to 
develcp those relationships or differences. FFT components 
9-14 have otvious differences; however, FFT components 20, 
24 and 25 also have large F-statistics but do not fall in 


the rejection area. These components may also be 
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TABLE 3 
F-STATISTICS OF FFT COMPONENTS 


i 

| 

| 

| 

| FFT AMONG WITHIN TOTAL 

| CCMPCNENT GROUP SUM GROUP SUM OF SQUARES F-STAT 
| 1 0.121 17.69 17.81 0.727 
2 0.112 18.48 18.59 0.642 
3 0.108 20.65 20.76 0.557 
4 0.135 25.40 25.54 0.562 
5 02 17> 32.74 32.91 0.558 
6 0.254 41.72 41.98 0.647 
7 0.422 50.10 50.52 0.893 
8 0.962 82.83 83.79 1.231 
9 7.591 390.2 397.8 2.063 
10 49.09 1681.0 1730.1 3.097 
11 53.01 173224 1785.1 3.246 
12 61.56 2079.7 2133.3 3.152 
13 45.73 1832.8 1878.6 2.646 
| 14 3.864 250.5 254.4 1.636 
| 15 0.539 98.13 98.67 0.582 
16 0.366 110.1 110.5 0.353 
17 0.504 151.0 151.5 0.354 
18 0.943 140.3 141.2 0.7.13 
| 19 1.299 158.5 159.8 0.370 
20 2.671 200.7 203.4 1.412 
21 1.853 226.3 228.2 0.868 
22 26277 236u1 234.3 1.041 
23 19842 197.6 199.3 0.919 
24 2.512 181.9 184.4 1.476 
25 2.144 171.3 173.5 1.327 
26 1.257 143.7 145.0 0.928 
| 27 1.233 113.9 115.1 1.149 
es 


an cs cs aS A ES A es Sa a eS A SE a A Se 


\ ees T= ae an mae | 
{ | 
{ FABLE. 3 \ 
| F-STATISTICS OF FFT COMPONENTS (con't) | 
| 28 0.858 VAS3'3 11607 S785 | 
! 29 0.547 G3 <.02 64.27 0.911 
| 30 0.356 5S als 56. 13 0.677 | 
31 0.183 41.41 By Wes) 0.469 
| SZ 0.051 13579 13.84 02389 
| S16 0.031 10... 92 10.96 On 94 
Ss £ ae ees ae 
c Te oS pe oe, a, aaa | 
| 
TABLE 4 i 
| F-STATISTICS OF THE MEAN VARIANCES, | 
| MEAN SKEWNESS AND MEAN KUJRTOSIS | 
| | 
| AMONG Wi THIN TOA { 
| GROUP SUM GROUP SUM OF SQUARSS F-STAT | 
VARIANCE 3. age ES Se LoD “Eb 4.070 E6 Dass 10 | 
SKEWNESS Seeed-E10c 703 -E41 WeoO2eintd r#2.2588 { 
| KURTOSIS 62250. BAT taeow (E13 Teo ets 2D a 0 | 
ee { 

eee ee a ee | 
interesting to examine. In addition, it should be noted 


tmiat ther statistics: of ..the FFT associated with each 
respective modulation technique also display some striking 
differences. These statistics may prove in some way to bea 


fingerprint of the modulation technigues themselves. 
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APPENDIX E 

RESULTS OF F-TEST 
F-STATISTICS FOR THE 1 COMPONENT 
SUM OF SUMS = 0.2091487882294994D 02 
SUM OF SUMS**2 = 0.1810655469453404D 02 
SUM OF SQUARE OF SUMS = 0.4128097510679048D 02 
AMCNG GROUP SUM OF SQUARES = 0.1211883136154515D 00 
WITHIN GROUP SUM OF SQUARES = 0.1769374494346614D 02 
TOTAL SUM OF SQUARES = 0.1781493325708159D 02 
DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.148500000C000000D O04 


F-STATISTIC = 0.7265063214387600D 00 


F-STATISTICS* FOR THE 2 COMPONENT 

SUM OF SUMS = 0.2305792717805020D 02 

SUM OF SUMS*¥*2 = 0.1894701234015291D 02 

SUM OF SQUARE OF SUMS = 0.4663825126655468D 02 

AMCNG GROUP SUM OF SQUARES = 0.1119371755000360D 00 

WITHIN GROUP SUM OF SQUARES = 0. 18480629827483736D 02 

TOTAL SUM OF SQUARES = 0.1859256700298739D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.21485000000000000D 04 
F-STATISTIC = 0.64247464050 14737D 00 


F-STATISTICS FOR THE 3 COMPONENT 

SUM OF SUMS = 0.2693734535827635D 02 

SUM OF SUMS**2 = 0.2124674921974186D 02 

SUM OF SQUARE OF SUMS = 0.5921566381358867D 02 

AMCNG GROUP SUM OF SQUARES = 0. 1084095881685183D 00 

WITHIN GROUP SUM OF SQUARES = 0.2065459258160598D 02 

TOTAL SUM OF SQUARES = 0.2076300216977449D 02 

DEGREES OF FREEDOM IN NUMERATOE = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.5567362242775778D 00 
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P-STATISTICS FOR THE 4 COMPONENT 

SUM OF SUMS = 0.3175688631627065D 02 

SUM OF SUMS**2 = 0.2621104542300214D 02 

SUM OF SQUARE OF SUMS = 0.8070513093726266D 02 

AMONG GROUP SUM OF SQUARES = 0. 1347180903696015D 00 

WITHIN GROUP SUM OF SQUARES = 0.2540399411362951D 02 

TOTAL SUM OF SQUARES = 0.2553871220399911D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 


F-STATISTIC = 0.5624997485041869D 00 

P=STATISTICS FOR THE 5 COMPONENT 

SUM OF SUMS = 0.3669104349531116D 02 

SUM CF SUMS**2 = 0.3380633740317695D 02 

SUM OF SQUARE OF SUMS = 0.1069682473356158D 03 

AMCNG GROUP SUM OF SQUARES = 0.1721940248396148D 00 
AITHIN GROUP SUM OF SQUARES = 0.327/3665492982079D 02 
TOTAL SUM OF SQUARES = 0.3290884895466041D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.14000090000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D O04 
F-STATISTIC = 0.5579331866788867D 00 


F—-STATIS£ZICS FOR THE 6 COMPONENT 

SUM OF SUMS = 0.4160097878948431D 02 

SUM OF SUMS**Z = 0.431323€642020094D 02 

SUM CF SQUARE OF SUMS = 0.1408215358007360D 03 

AMCNG GROUP SUM OF SQUARES = 0.2544544005119440D 00 

WITHIN GROUP SUM OF SQUARES = 0.4172415106219358D 02 

TCIAL SUM OF SQUARES = 0.4197860546270552D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 92 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.6468757561623442D 00 
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F-STATISTICS FOR THE 7 COMPONENT 

SUM CF SUMS = 0.4810387092348111D 02 

SUM OF SUMS**2 = 0.5206409381611744D 02 

SUM OF SQUARE OF SUMS = 0.1964653659903027D 03 

AMONG GROUP SUM OF SQUARES = 0.4219987280210724D 00 

WITHIN GROUP SUM OF SQUARES = 0.5009944015621441D 02 

TOTAL SUM OF SQUARES = 0.5052143888423548D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.8934632362547178D 00 


F=STATISTICS FORALAE 8 COMPONENT 

SUM OF SUMS = 0.5623255941290002D 02 

SUM CF SUMS**2 = 0.8590285622098203D 02 

SUM CF SQUARE OF SUMS = 0.3070018365343552D 03 

AMCNG GROUP SUM OF SQUARES = 0.9619512065933320D 00 

WITHIN GROUP SUM OF SQUARES = 0.8283283785563848D 02 

TOTAL SUM OF SQUARES = 0.8379478906223181D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.1231824736914027D 01 


F-STATISTICS FOR THE 9 COMPONENT 
SUM OF SUMS = 0.9931605295569187D 02 

SUM OF SUMS#*2 = 0.4043531056175881D 03 

SUM OF SQUARE OF SUMS = 0.1416642417283150D 04 

AMCNG GROUP SUM OF SQUARES = 0.7590638589699638D 01 

WITHIN GROUP SUM OF SQUARES = 0.3901866814447566D 03 

TOTAL SUM OF SQUARES = 0.3977773200344562D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000006000000D 04 
F-STATISTIC = 0.2063499133280512D 01 
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F-STATISTICS FOR THE 10 COMPONENT 

SUM OF SUMS = 0.2029942628781924D 03 

SUM OF SUMS¥*¥*2 = 0.1757564757819924D 04 

SUM OF SQUARE OF SUMS = 0.7655692250451960D 04 

AMCNG GROUP SUM OF SQUARES = 0.4908580866354581D 02 

WITHIN GROUP SUM OF SQUARES = 0.1681007835315404D 04 

TOTAL SUM OF SQUARES = 0.1730093643978950D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.3067309684192641D 01 


PSSTATESTICS ‘FORSTHE 11 COMPONENT 

SUM OF SUMS = 0. 2004773302247246D 03 

SUM CF SUMS**2 = 0.1811875077561617D O04 

SUM OF SQUARE OF SUMS = 0.7980439863108266D 04 

AMCNG GROUP SUM OF SQUARES = 0.5301029200839379D 02 

WITHIN GROUP SUM OF SQUARES = 0.1732070678930534D 04 

TOTAL SUM OF SQUARES = 0.1785080970938928D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
PoOTAPESTLC = <053.2:4633 2537532912) (01 


F-STATISTICS FOR THE 12 COMPONENT 

SUM OF SUMS = 0.£2210834652773109D 03 

SUM CF SUMS¥**2 = 0.2165859765218944D 04 

SUM OF SQUARE OF SUMS = 0.9414279932820438D 04 

AMCNG GROUP SUM OF SQUARES = 0.6155753358218842D 02 

WITHIN GROUP SUM OF SQUARES = 0.2071716965890740D 04 

TOTAL SUM OF SQUARES = 0.213327449 9472928D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.3151731454585565D 01 
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F-STATISTICS FOR THE 13 COMPONENT 

SUM OF SUMS = 0.1979583434471791D 03 

SUM OF SUMS**2 = 0.1904701075347160D 04 

SUM OF SQUARE OF SUMS = 0.7185446672840862D 04 

AMCNG GROUP SUM OF SQUARES = 0.4572946290150773D 02 

WITHIN GROUP SUM OF SQUARES = 0. 1832846608618751D 04 

TOTAL SUM OF SQUARES = 0.1878576071520259D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D O04 
F-STATISTIC = 0.26464786714598604D 01 


F-STATISTICS FOR THE 14 COMPONENT 

SUM CF SUMS = 0292277362726 14619D 02 

SUM OF SUMS**2 = 0.2600827196886254D 03 

SUM OF SQUARE OF SUMS = 0.9540726819081072D 03 

AMCNG GROUP SUM OF SQUARES = 0.3863985704619236D 01 

WITHIN GROUP SUM OF SQUARES = 0.250541992869 5444D 03 

TOTAL SUM OF SQUARES = 0.2544059785741636D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.1635887377498235D 01 


F=STATLSTICS, FORTHE 45, COMPONENT 

SUM OF SUMS = 0.6880926336210145D 02 

SUM CF SUMS¥**2 = 0.1018300362534556D 03 

SUM OF SQUARE OF SUMS = 0.3695062863415098D 03 

AMCNG GROUP SUM OF SQUARES = 0.5385863804584070D 00 

WITHIN GkOUP SUM OF SQUARES = 0.9813497339004G56D 02 

TOTAL SUM OF SQUARES = 0.9867355977049896D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1400000000000000bD 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.14850000000000900D O04 
F-STATISIIC = 0.58219339 710755960, 00 
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F-STATISTICS FOR THE 16 COMPONENT 
SUM OF SUMS = 0.6961387326450431D 02 

SUM OF SUMS**2 = 0.1137157033969463D 03 

SUM OF SQUARE OF SUMS = 0.3597053597878584D 03 

AMCNG GROUP SUM OF SQUARES = 0.3663260306209388D 00 

WITHIN GROUP SUM OF SQUARES = 0.1101186497990677D 03 

TOTAL SUM OF SQUARES = 0.1104849758296886D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.3528623485827816D 00 


F—STATISTICS FOR THE 17 COMPONENT 

SUM OF SUMS = 0.7052854165691613D 02 

SUM OF SUMS**2 = 0.1548042004914454D 03 

SUM OF SQUARE OF SUMS = 0.3820324148469269D 03 

AMONG GROUP SUM OF SQUARES = 0.5041406896350342D 00 
WITHIN GROUP SUM OF SQUARES = 0.1509838763429761D 03 
TOTAL SUM OF SQUARES = 0.1514880170326111D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D O4 
F-STATISTIC = 0.3541763825767676D 00 


F-STATISTICS FOR THE 18 COMPONENT 

SUM OF SUMS = 0.7950452893326030D 02 

SUM CF SUMS**2 = 0. 1454127656122291D 03 

SUM OF SQUARE OF SUMS = 0.5157161626920556D 03 

AMCNG GROUP SUM OF SQUARES = 0.9431815463208067D 00 

WITHIN GROUP SUM OF SQUARES = 0. 14025560398539086D 03 

TOTAL SUM OF SQUARES = 0.14119878553 16294D 03 

DEGREES OF FREEDOM IN NUMEKATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.21485000000000000D 04 
F-STATLSTIC = 0.7133020797581566D 00 
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F-STATISTICS FOR THE 19 COMPONENT 

SUM OF SUMS = 0.8206200282896035D 02 

SUM OF SUMS**2 = 0.1642503769158388D 03 

SUM OF SQUARE OF SUMS = 0.5788746875066006D 03 

AMCNG GROUP SUM OF SQUARES = 0.1299298669532475D 01 

WITHIN GROUP SUM OF SQUARES = 0. 1584616300407728D 03 

TOTAL SUM OF SQUARES = 0.1597609287103052D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMZNATOR = 0.1485000000000000D 04 
F=STATISTIC = 0586972768097 11277). 00 


F-STATISTICS FOR THE 20 COMPONENT 

SUM OF SUMS = 0.9520048302882186D 02 

SUM OF SUMS**2 = 0.2094086558414699D 03 

SUM OF SQUARE OF SUMS = 0.8713219218815827D 03 

AMCNG GROUP SUM OF SQUARES = 0.2671131239535162D 01 
WITHIN GROUP SUM OF SQUARES = 0.2006954366226541D 03 
TOTAL SUM OF SQUARES = 0. 2033665678621893D 03 

DEGREES OF FREEDOH IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.1411744637781583D 01 
F-STATISTICS FOR THE 21 COMPONENT 


SUM OF SUMS = 0.8994139952150146D 02 

SUM OF SUMS¥**2 = 0.2335862791765991D 03 

SUM OF SQUARE OF SUMS = 0.7245755635510632D 03 

AMCNG GROUP SUM OF SQUARES = 0.1852785403586404D 01 

WITHIN GROUP SUM OF SQUARES = 0.2263405235410884D 03 

TOTAL SUM OF SQUARES = 0.2281933089446748D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.148500000C000000L 04 
F-STATISTIC = 0. 8682828488687504D 00 
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F-STATISTICS FOR THE 22 COMPONENT 

SUM OF SUMS = 0.1035648300735959D 03 

SUM OF SUMS**2 = 0.2414810868436007D 03 

SUM OF SQUARE OF SUMS = 0.9427341493765320D 03 

AHCNG GROUP SUM OF SQUARES = 0.2276892141650119D 01 

WITHIN GROUP SUM OF SQUARES = 0.2320537453498354D 03 

TOTAL SUM OF SQUARES = 0.2343306374914855D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.1040764077320930D 01 


P=STATISTICS FOR THE 23 COMPONENT 

SUM OF SUMS = 0.9031234717154746D 02 

SUM OF SUMS**2 = 0.2047768116906234D 03 

SUM OF SQUARE OF SUMS = 0.7149379755971323D 03 

AMCNG GROUP SUM OF SQUARES = 0.1711833054881912D 01 

WITHIN GROUP SUM OF SQUARES = 0.1976274319346€521D 03 

TOTAL SUM OF SQUARES = 0.1993392649895340D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02. 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F=STATISTIC = 0.9 1878 22552243536) 00 


F-STATISTICS FOR THE 24 COMPONENT 

SUB OF SUMS = 0.9171706807335522D 92 

SUM OF SUMS**2 = 0.1900438745495707D 03 

SUM CF SQUARE OF SUMS = 0.8139562480340389D 03 

AMCNG GROUP SUM OF SQUARES = 0. 2531548763025405D 01 

WITHZN GROUP SUM OF SQUARES = 0.1819043120692303D 03 

TOTAL SUM OF SQUARES = 0. 1844358008322557D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D JZ 
DEGREES CF FREEDOM IN DENOMINATOR = 0.148500000000000uUD 04 
F-STATISTIC = 0.1476188171339999D 01 
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F-STATISTICS FOR THE 25 COMPONENT 

SUM OF SUMS = 0.8440126623444366D 02 

SUM OF SUMS**2 = 0. 1782254895437147D 03 

SUM OF SQUARE OF SUMS = 0.6892901073945840D 03 

AMCNG GROUP SUM OF SQUARES = 0.2143851912627546D 01 

WITHIN GROUP SUM OF SQUARES = 0.1713325884697689D 03 

TOTAL SUM OF SQUARES = 0.1734764403823964D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.1327251499840134D 01 


F-STATISTICS FOR THE 26 COMPONENT 

SUM OF SUMS = 0.6923579472492147D 02 

SUM OF SUMS**2 = 0.1481985058366125D 03 

SUM OF SQUARE OF SUMS = 0.4452980334 185913D 03 

AMCNG GROUP SUM OF SQUARES = 0.1257250153391605D 01 

WITHIN GROUP SUM OF SQUARES = 0.1437455255024266D 03 

TOTAL SUM OF SQUARES = 0.1450027756558182D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGKEES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.9277389287476910D 00 


F-STATISTICS FOR THE 27 COMPONENT 

SUM OF SUMS = 0.6366092822129744D 02 

SUM OF SUMS**2 = 0.1177906819912730D 03 

SUM OF SQUARE OF SUMS = 0.3935159041535181D 03 

AMCNG GROUP SUM OF SQUARES = 0.1233349853537058D 01 

WITHIN GROUP SUM OF SQUARES = 0. 1138555229497378D 03 

TOTAL SUM OF SQUARES = 0.1150888728032749D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.1149027974258136D 01 
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F-STATISTICS FOR THE 28 COMPONENT 

SUM OF SUMS = 0.5387592494868525D 02 

SUM OF SUMS**2 = 0.1186133442790129D 03 

SUM OF SQUARE OF SUMS = 0.2792595957915495D 03 

AMCNG GROUP SUM OF SQUARES = 0.8575190985312517D 00 

WITHIN GROUP SUM OF SQUARES = 0.1158207483210974D 03 

TOTAL SUM OF SQUARES = 0.1166782674196286D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.7853366268738306D 00 


F-STATISTICS FOR THE 29 COMPONENT 

SUM OF SUMS = 0.4536144151906700D 02 

SUM OF SUMS**2 = 0.6564305139058856D 02 

SUM OF SQUARE OF SUMS = 0.1919053830677748D 03 

AMONG GROUP SUM OF SQUARES = 0.5472802462192574D 00 

WITHIN GROUP SUM OF SQUARES = 0.6372399755991081D 02 

TOTAL SUM OF SQUARES = 0.6427127780613007D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.9 109723144851782D 00 


F-STATISTICS FOR THE 30 COMPONENT 
SUM OF SUMS = 0.3732568029970994D 02 

SUM OF SUMS**2 = 0.5705765228634419D 02 

SUM OF SQUARE OF SUMS = 0.1284820507967475D 03 

AMCNG GROUP SUM OF SQUARES = 0.3560162347433721D 00 

WITHIN GROUP SUM OF SQUARES = 0.5577283177837672D 02 

TOTAL SUM OF SQUARES = 0.5612884801312009D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.14000000000000G0D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.6770886363437515D 00 
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F-STATISTICS FOR THE 31 COMPONENT 

SUM OF SUMS = 0.2973876286269883D 02 

SUM OF SUMS**2 = 0.4217833610659068D 02 

SUM OF SQUARE OF SUMS = 0.7726205204918280D 02 

AMONG GROUP SUM OF SQUARES = 0.18302 45094226048D 00 

WITHIN GROUP SUM OF SQUARES = 0.4140571558609885D 02 

TOTAL SUM OF SQUARES = 0.4158874009552145D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.4688645251806333D 00 


F-STATISTICS FOR THE 32 COMPONENT 
SUM OF SUMS = 0.2142367689889270D 02 

SUM OF SUMS¥**2 = 0.1414820425416038D 02 

SUM OF SQUARE OF SUMS = 0.3565723402182890D 02 

AMONG GROUP SUM OF SQUARES = 0.5058971897285643D-01 

WITHIN GROUP SUM OF SQUARES = 0.1379 163191394209D 02 

TCTIAL SUM OF SQUARES = 0.13842221632914S55D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 62 
DEGREES OF FREEDOM IN DENOMINATOR = 0.14850000000000005 04 
F-STATISTIC = 0.38908548284%65455D 00 


F-STATISTICS FOR THE 33 COMPONENT 
SUM OF SUMS = 0.1732004830070321D 02 

SUM OF SUMS**2 = 0.1115582097066481D 02 

SUM OF SQUARE OF SUMS = 0.2312899538600082D 02 

AMCNG GROUP SUM OF SQUARES = 0.3130057176754676D-01 

AIIHIN GROUP SUM OF SQUARES = 0.1092453101680480D 02 

TOTAL SUM OF SQUARES = 0. 1095583158857235D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.140000C000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.3039120267386334D 00 
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F-STATISTICS OF MEAN VARIANCE, SKEWNESS AND KURTOSIS 


F-STATISTICS FOR THE 1 COMPONENT 

SUM CF SUMS = 0.1111617439422250D 05 

SUM OF SUMS**2 = 0.4319425777749392D 07 

SUM CF SQUARE OF SUMS = 0.1861260733960752D 08 

AMONG GROUP SUM OF SQUARES = 0.3143833877400904D 06 

WITHIN GROUP SUM OF SQUARES = 0.3755407373518861D 07 

TOTAL SUM OF SQUARES = 0.4069790761258952D 07 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.48000000000000000 05 
F-STATISTIC = 0.2870223636519026D 01 


R-STATISTiICS EOR THE 2 COMPONENT 

SUM OF SUMS =-0. 3062348822555926D 07 

SUM CF SUMS**2 = 0.7751235739069674D 12 

SUM OF SQUARE OF SUMS = 0.2376364181020636D 13 

AMONG GROUP SUM OF SQUARES = 0.5306562101878761D 11 

WITHIN GROUP SUM OF SQUARES = 0.7031125381184634D 12 

TOTAL SUM OF SQUARES = 0.7561781591372509D 12 

DEGREES OF FREEDOM IN NUMERATOR = 0.14G60000000000000D 02 
DEGREES CF FREEDOM IN DENOMINATOR = 0.4800000000000000D 03 
F-STATISTIC = 0.2587626620217669D 01 


P-STATISTICS FOR THE 3 COMPONENT 

SUM OF SUMS = 0.1087692759375071D 08 

SUM OF SUMS**2 = 0.1368235571001090D 14 

SUM CF SQUARE OF SUMS = 0.3521258942662316D 14 

AMCNG GROUP SUM OF SQUARES = 0.8280430048704694D 12 

WITHIN GEOUP SUM OF SQUARES = 0.1261530754556778D 14 

TOTAL SUM OF SQUARES = 0. 1344335055043825D 14 

DEGREES GF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.4800000000000000D 03 
F-STATISTIC = 0.2250444214596077D 01 
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